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ABSTRACT 



This research is a study in seepage analysis for engineered 
structures built below the groundwater table; in particular an 
analysis of under drain systems for such structures. General 
seepage theories and design considerations are discussed as 
part of the literature review. Also design concepts and 
methods of analysis are covered which gives a brief overview of 
design considerations and techniques which are currently used 
•for this type of seepage problem. 

Key variables to the design of an underdrain system are studied 
in detail by using a two dimensional, steady state finite 
element model. Such variables studied are drain spacing, drain 
width, anisotropic hydraulic conduct i vi ty , thickness of the 
aquifer, head differential and sloping aquifers. 

Results are presented in the form of dimesionless plots and 
nomographs for drain discharge and free surface plots for the 



different variables studied. 
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CHAPTER ONE 



INTRODUCTION 

Far civil engineered structures which interact with soil, a 
major consideration -for design is the presence o-f groundwater. 
Dams, buildings, dry dock -facilities, underground structures, 
pavements and earth retaining walls must be designed 
considering the quantity of seepage, the position of the 
free surface, seepage forces and seepage control. Good 
drainage incorporated in the design of these and other 
structures can often benefit through substantial cost 
reductions as well as increasing the safety and usefulness of 
the structure. 



Many embarassing and sometimes disastrous problems often result 
from man's endeavors to control seepage. Saturated soil 
backfills have caused many retaining wall failures due to 
increased hydrostatic pressures following extensive rainfall. 
Highways and roadways prematurely deteriorate and fail due to 
trapped water in the pavement section. Basement slabs, pond 
liners and drydock facilities have been severely damaged by 
hydrostatic uplift pressures and seepage. Most dam and 
slope failures can be attributed to excessive seepage or 
uplift pressures. These examples and untold others indicate 
that good seepage and drainage control together with an adequate 



understanding of groundwater is a key element in the design of 
most civil engineering structures. 

Many civil engineering structures are located in areas where 
groundwater causes problems of uplift and seepage. The design 
therfore must consider all problems associated with high 
groundwater and what measures must be ~ taken to assure the 
integrity of the overall design. Casagrande ( 1 965) , in talking 
about engineering risks, concludes that, "The margin of safety 
that we incorporate into our structures should b^ar a direct 
relationship to the magnitude of potential losses, and it must 
also take into account the range of uncertainty involved." 

This important concept hopefully will be considered throughout 
the presentation of facts and data obtained in this paper. 

Prior to the 20th century, analysis of seepage and the control 
of drainage for enginerring structures was done by experienced 
engineers who based their designs on trial and error methods 
which sometimes were successful and sometimes not. Though 
experience has proven to be very neccesary and helpful in the 
design of drainage systems, very useful tools employed today 
are numerical methods such as finite element and finite 
difference techniques which require the use of the digital 
computer. Their problem solving capabilities have made 
possible solving very complicated and complex groundwater 
pr obi ems. 
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The purpose of this research is to study the design of 
underdrain systems -for large structures placed below the 
groundwater table. In chapter two, basic seepage theories and 
design requirements, designs concepts and useful methods of 
analysis are presented as a brief literature review. This 
chapter covers critical aspects of designs that should be 
considered and contains a brief review of methods used by 
geotechnical engineers. 

In chapter three, the finite element program developed by Dr. 
M. Aral, of the School of Civil Engineering at Georgia 
Institute of Technology, is discussed describing how the 
design of building underdrain systems can be modeled using 
this technique. This finite element program made possible 
studying effects of key variables such as permeability, 
hydraulic head, drain geometry, depth of groundwater, and 
anisotropic soils. Program verification is performed using 
another program, SEEP84, modified by the Civil Engineering 
Department at Virginia Polytechnic Institute and State 
University, for the microcomputer. Dupuit’s solution is also 
used to compare selected cases with the finite element 
pr ogr ams. 
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Chapter four presents data generated from the finite element 
study such as position of the free surface and flow quantity 
for selected systems of underdrains suitable for use beneath a 
foundation- The data is presented in the form of charts and 
nomographs to be used as a guide for underdrain system design. 

Chapter five summarises the results of the study and presents 
appropriate conclussions. 
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CHAPTER 2 



LITERATURE REVIEW 



2.1 Seepage Theories and Design Requirements 



G r ■ o u nd water 



Soils, sometimes described as a porous media, is a mass made up 
of solid matter having interconnected voids, passages, fissures 
or pore spaces. These pore spaces can be filled by a liquid 
and or a gaseous phase. Figure 2.1 shows several typical soil 
and rock indices which depicts how pore spaces can act as 
conduits for water or subcapillary openings in which water may 
be held by adhesive capillary forces (Meinzer, 1942) , 

The soil mass forms a complex network of irregular passages 
through which the fluid may flow. A key parameter of 
groundwater movement is the porosity of the soil. As gravity 
acts on the fluid, it moves through the pore spaces at varying 
rates through each open pore space. From a macro scale, this 
flow has been mathematically described by several authors 
(Bear, 1981; Terzaghi, 1943; Cedergren, 1977; Harr, 1962; 

Bouwer, 1973; Freeze, 1971). 
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Figure 2.1 Diagram showing soil and rock interstices. A., Well 
sorted sand deposit having high porosity; B. Poorly sorted 
deposit having low porosity; C. Particles that maybe porous 
themselves; D. Deposit whose porosity may be diminished by 
deposition of mineral matter in the interstices; E. Rock that 
is porous by solution; F. Rock that is porous by fracturing 
(after Meins er, 1942) 
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Moisture content in sails can vary with depth -from the surface, 
as shown by Figure 2„2. Most seepage problems considers the 
zone of saturation which is the zone where the soil mass is 
completely saturated,. The terms groundwater level, phreatic 
water level, water table, free surface or groundwater surface 
are all synonymous to denote the level of water where 
atmospheric pressure exists. 



( iround surface 
Soil water /one 

A 

•' 1 77 ” ' r 




Soil w.iter 



Well 



Intermediate 
(Vadose water) 
/one 



T 




Pellicular and 
gravitational water 

litli!fiillll!il!iLv.Jlilllllll!lim 

Capillary water 

Capillary fringe Groundwater 
^ Impervious 




Water- 
table 






Figure 2.2 



Subsurface groundwater distribution, 
(after Bear, 1979) 
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In dealing with capillary fringe in fine grained soils the 
height of water can very significant. As the pore diameter 
decreases, the capillary rise increases. Although the soil i 
considered unsaturated, the soils moisture content will be 
higher with depth. 

The zone of capillary rise and fringe zone may contain large 
quantities of water for fine grained soils. Water removal by 
drainage, due to capillary forces, may be difficult for such 
fine grained soils. Figure 2.3 shows how the actual moisture 
distribution is appr ox i mated by a step distribution. This 
approximation of the capillary fringe level above the free 
surface, h c , is often neglected in some seepage problems. 
Several emperical methods, based on grain size diameter and 
porosity of soil, are available to estimate h c (Mavis and 
Tsui, 1939; Pol ubari nova, 1952, 1962; Silin Bekchurin, 1958). 
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Figure 2.3 Approx i mat i on of phreatic surface and capillary 
fringe zone. (after Bear, 1979) 



Many textbooks on seepage and groundwater discuss water bearing 
and retention capibilities of the capillary fringe and 
movement of moisture through this zone (Terzaghi, 1943; Kovacs, 
1981; Sowers, 1979). However, when this height of water is ver 
small relative to the total thickness of the aquifer, the 
capillary fringe is sometimes neglected. The geotechnical 
engineer should be concerned with the unsat urated-saturated 
flow condition when the stability of the engineered structure 
is to be analyzed, such as a dam. 



Groundwater analysis for seepage problems must consider the 
quantity of water to be drained, the soils resistance to 
drainage, and what effects the drainage may have on the soil 
and the structure. Movement of groundwater can be categorized 
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either as seepage; capillary rise or flow; percolation; or 
turbulent flow. For this study, only movement of groundwater 
by seepage is considered. 

Design Requirements 

Structures built below the watertable must be designed 
considering the following: 

(1) structural ability to withstand the applied hydrostatic 
water pressures, 

<2) drainage control which short-cuts the flow or movement 
of groundwater from getting to the structure, 

(3) drainage control which is incorporated with the overall 
design to remove the groundwater, 

(4) assure water velocities be small enough to prevent 
excess particle migration and erosion. 

Savings of dollars, time and natural resources can be 
recognised when structures can be considered on dry, stable 
material;in addition to thinner basement slab thicknesses may 
be used. 



10 



Darcy" s Law 



In 1856, Henry Darcy investigated the -flow of water in vertical 
homogeneous sand. This work has developed the most fundamental 
relationship which is still used in seepage analysis today. 
Darcy's Law states: 



Q = k i A (2.1) 

where Q <lVt) is a volume discharge or flow of a given 

2 

cross-sectional area, A (L ) in a given time. The discharge is 
also a function of the hydraulic gradient, i (L/L) and a 
resistance coefficient, k (L/T) , the permeability of the soil. 
Investi gators, Muskat (1937), Taylor (1949) and Leonards 
(1962), have concluded that Darcy's Law is valid where 
velocities are low or when flow is laminar, i.e. Reynolds 
number less than 10. 

Hydraulic Conductivity 



The soils resistance to flow or the coefficient of 
porporti onal i ty , k, for Darcy’s Law, is also called the 
hydraulic conduct i vi ty . Hydraulic conductivity is a 
coefficient which expresses the ease or resistance the water 
has flowing through the soil. It depends on the soils matrix 
and fluid properties of the permeant. 
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Many researchers have performed tests to attempt, to develops 
empirical relationships of the soils properties to hydraulic 
conductivity. Borne of these properties considered are grain or 
pore size distribution, shape of the grains or pores, 
tortuosity, specific surface and porosity (Mitchell, 1976$ Kozeny, 
1927; Carman, 1956; Cedergren, 1977; Kovacs, 1981). All report 
that careful field i nvesti agat i on and careful conductivity 
studies are necessary for each design. Figure 2,. 4 is presented 
to show some characteristics of various soil types and most 
important, methods for determination of selected properties 
(Terzaghi and Peck, 1948). 

In the unsaturated zone, the hydraulic conduct i vi ty , moisture 
content and pressure head is extremely complex and not well 
understood. Again, field measurement is the best method for 
obtaining the best parameters, such as the use of theistor 
psychrometer (Kay and Low, 1970) or gamma ray absorption 
techniques (Ferguson, 1970; Aral, 1983). However, many 
empirical formulas have been developed to relate the soil 
parameters to a relative hydraulic conductivity of the 
un-saturated zone (Kovacs, 1981; Irmay, 1954; Childs and 
Col 1 i s-George, 1950; Gardner, 1958; Brooks and Corey, 1966; 

Maslia, 1980; Reeves and Duriguid, 1975). 
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Coefficient of Permeability k in cm per sec (log scale) 
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Figure 2.4 Permeability and Drainage Characteristics for Various Soil Types 



2.2 Design Concepts 



Seepage Forces 

The effects of water on the soil and the structure is 
controlled by the pressures in the water. In soil mechanics, 
the water pressure is defined as neutral stress. For water not 
moving, or steady state, the pressure can be computed from the 
basic equation of hydrostatics 



h = 




( 2 . 2 ) 



where h, energy head, is equal to the position, z, plus the 
pressure head or piezometric head, F'/^ . Typically in 

seepage designs, the velocity head is very small and may be 
neglected. Then the pressure head plus elevation is equal to 

head, h (Cedergren, 1977). 

For the condition that water pressure is greater than the total 
load, say of the soil and structure, the pressures are termed 
uplift pressures and liquifaction can occur. Therfore the 
control of uplift pressures must be a primary design 
consideration. 
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Critical Velocity 



According to Darcy’s Law, the discharge velocity of seepage 
through soil is 



v , = k x i 

a 



( 



) 



Seepage is only occur! ng through the pores and voids. 

Therefore the true seepage velocity should be a function of the 
porosity of the soil mass. Re-writing equation (2.3) the 
seepage velocity can be expressed as 



v 



sn 



k x i 




( 2 . 4 ) 



where i\ e is the effective porosity. 

In flow conditions, the saturation of soil is moving into or 
away from the unsaturated soil. The capillary head, h c , is 
working with gravity to increase the hydraulic gradient. For 
this condition the hydraulic gradient can be expressed 

h + h 

i = , c 

L (2.5) 

Figure 2.5 depicts this condition for flow through a sand. 
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Reservoir 





(b) 



Figure 2.5 Force of saturation. 



Boiling occurs when seepage pressures in an upward direction 
exceeds the downward -force o-f the soil. This must be analysed 
for all seepage designs. Critical hydraulic gradient, , can 

be expressed 



where 




i c 






( 2 . 6 ) 



total unit weight of soil 
unit weight of water- 

buoyant unit weight of so.i 1 (NAVFAC, DM-7. 1 ) . 
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Fi 1 ter s 



A very important design consideration is the possibility of 
seepage -forces moving erodible soils and rock, termed piping. 
Betram (1940), under Terzaghi and Casagrande, established the 
following criteria for filter design,, bashed on grain size. 



D (of filter) D (of filter) 

— ^ < 4 to 5 < — 

D (of soil) D (of soil) (2.7) 

oo 15 



This relationship is represented by Figure 2.6 to illustrate 
the prevention of piping by filters (Cedergren, 1977). 



LEGEND 




Figure 2.6 Illustration of prevention of piping by filters. 

(after Rankilor, 1981) 
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The ratio of D^5 of a filter to the Dg^ of a soil is 
considered the piping ratio. Many authors and organ i sat i ons 
have established other criteria for filters, similar to that of 
equation (2.7) (U.S. Army Corps of Engineers, 1941, 1955; U.S 

Bureau of Reclamation, 1955; Sherard, etal , 1963; Cedergren, 

1977; Thani kachal am and Sakth i vad i vel , 1974). 

Filters must be designed to meet two major requirements! (1) 
they must be safe with respect to erosion and piping, and (2) 
must have sufficient discharge capi bi 1 1 i t i es to remove the 
seepage quickly and efficiently. This must be done without 
inducing high seepage forces or hydrostatic pressures. The 
right side of equation (2.7) states that the D size of a 
filter should be at least 4 to 5 times the D-^ of the soil. 

This insures that the hydraulic conductivity of the filter and 
drains is sufficient to prevent the build-up of large seepage 
forces and hydrostatic pressures. 

There is a wide range of engineering properties of materials to 
be considered for filter drainage systems designs. Properties 
which may influence hydraulic conductivity are grain size, 
density, mi ner al og i cal composition, the nature of the permeant or 
the degree of saturation. Much work has gone toward methods 
for determination and factors which influence hydraulic 
conductivity (Milligan, 1975; Lambe and Whitman, 1969; 

Yemington, 1963; Bouwer, 1978). There is still much research 
necessary to understand better this property of soil. 
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Trench Drains 



For the application of drainage -for structures built below the 
water table, the use of trench drains, which acts as collectors 
of the seepage or may contain collector pipes, is a very 
efficient method. Figure 2.7 shows the typical cross-sect i on 
which may be considered for the application of this type of 
drainage system. Trench drains will be discussed further in 
the report and their utilization for this application will be 
studied in more detail. 




Figure 2.7 Typical application of a trench drain system. 
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Geosynthet i cs 



Today, new synthetic membranes are used to serve the purpose o-f 
either -filtration or seperation. This two -factors are very 
important parameters for most seepage problems. Figure 2.8 
shows comparisons of some of the features of membrane-wrapped 
drains to conventional granular pipe drains. Also Figure 2.9 
shows the bridging effect caused by good seperation of the 
larger materials and smaller particles. 



Granular drain 



Membrane-wrapped drain 
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Pipe drains 
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rodding points. 
Expensive. 



Pipe must be 
levelled by 
engineer, to 
be an exact 

gradient. This is expen- 
sive. Trench must be left 
open for this, and can get 
silted up during storms. 

Well-point 
pumping 
may be 
necessary to 
lay the pipe. 
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Sand and stone must be 
supplied and laid to a strict 
specification. Very expensive. 

Pipe is expensive to buy, and 
must be laid by hand. Trench 
must therefore be wide, and 
timbered if deep. 
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The deep transporter means 
that settlement can take place 
without the drain blocking — 
unlike a pipe 



Wide trench 
means more 
excavation and 
more backfill than is 
technically necessary 



Can be 
taken deeper 
because sides 
do not need 
support 




t h 




Cost of shoring 
up open trench 
is eliminated 

Granular fill only 
supports trench, 
so can be cheap 
specification 

This drain is very 
much quicker to 
construct & can 
be laid in flooded 
ground. 

Narrower 

width because no pipe to be 
laid, so men do not have to 
go into the trench. 



Figure 2.8 Conventional granular pipe trench drain and 
membrane-wrapped trench drain. (after Rankilor, 1981) 
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1 1 





Magnified cross section of permeable membrane 

filter drain 




Filter Cake) 



Figure 2.9 Bridging e-f-fects caused by use o-f membrane -filter 
material. (after Eankilor, 1981) 
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A good representation of how geosynthet i cs combine action of 
separation and filtration, is shown by Figure 2.10. As water 
has more area of flow with the membrane, the seepage can be 
considered to be at a faster rate or drain quicker (Rankilor, 
1981) 



Without membrane 



With membrane 





O q o o “d 



y 

/ w 



Cui 



With membrane 



Large surface area 
for moisture loss 



t 



Small 
stone/clay 
contact area 



Figure 2.10 Seperation effects by using a membrane. 

(after Rankilor, 1981) 
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Methods of Analysis 



In analyzing seepage problems; geometry, boundary conditions, 
soil parameters, flow quantities and location of the free 
surface must be quantified by some means. Though this involves 
many consi derat i ons, the engineer must gain reasonable safe 
results for the design. 

Through the history of seepage analysis many techniques and 
methods have been developed, such as graphical, analytical and 
numerical solutions. Still used today is the graphical method 
of flow net construction. However with the complex geometries 
envoi ved and via the speed of computers, numerical models have 
been developed as very useful tools for the engineer. 

Flow Nets 



Many researchers have developed methods for analyzing flow 
through porous media. The oldest, developed by Casagrande and 
Terzaghi (1937), is the use of flow nets, a graphical approach 
to seepage analysis. 
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The continuity equation in the -form 



3u 


3v 


3w 






(2. 7) 


h 

3x 


3y 


+ — 

3z 


— 


0 





shows -for a Cartesian coordinate system, the quantity of water 
entering an element must equal that leaving. The terms u, v 
and w are discharge v-elocity components of directions x , y and 
2. respectfully. 



According to Darcy's Law, expressed by equation (2.3), the 
component of discharge velocity can be expressed 



u - -k 



3h 



3 x 



v = -k 



3 h 



3Y 



w 



ah 
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( 2 . 8 ) 



Substituting into equation (2.7) 



-k ( 3h/3x ) -k < 3h/3y ) -k(3h/3z) (2.9) 

a + a + a = o 



3 x 



9y 



a z 
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If k is constant for the element and then expressed in a 



typical two-dimensional form, equation (2.9) becomes the common 
Laplace equation. 



3 2 h 

2 + 

ax 



a 2 h 



ay 



= o 



( 2 . 10 ) 



The general conditon of -flow 
represented by Figure 2.11. 
presented in many textbooks 



analysis using -flow nets is 
Rules and good examples are 
(Cedergren, 1977; Harr, 1967). 



Free water surface 




Figure 2.11 General condition -for flow lines and 
equipotenti al lines in flow net construction. 
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Dupuit Assumptions 



For uncon-f ined seepage problems the engineer must calculate 
the location of the -free surface. Many problems of flow have 
been approximated based on the Dupuit theory of unconfined 
flow. Dupuit (1893) assumed; (1) for small inclinations of 
the line of seepage, the streamlines can be taken as horizontal 
and equal potential lines vertical, and (2) the hydraulic 
gradient was equal to the slope of the free surface and 
.invarient with depth (Harr, 1962). 

Harr (1962) gives many example problems and has derived the 
condition of free surface elevations as shown by Figure 2.12. 




Figure 2.12 Simple free surface notation. 
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In the case of Figure 2.12 the Dupuit formula is expressed as 



q 



k 




( 2 . 11 ) 



Many other reports and textbooks have presented charts and 
nomagraphs based on these assumptions. Later in this study a 
nomograph for a simple trench drain with blanket drain 
solution developed by the Federal Highway Admi n i strat i on (1980) 
will be presented. 



Numerical Analysis 

Flow through porous media has been approximated by many authors 
using ordinary non-linear differential equations and partial 
differential equations. Coupled with the speed of the computer, 
numerical solutions to many seepage problems have been 
developed. It should be pointed out that complex mathematical 
problems solved by numerical methods, are only as accurate as 
the math model can represent the conceptions and phenomina of 
seepage flow. In most cases, simplification is the rule, and 
numerical models can only be treated as approximations to the 
actual occurances. Nevertheless, the number crunching 
capibilities and time savings involved gives the engineer much 
needed time for the true analysis of the seepage problem. 

The analytical and graphical methods mentioned previously are 
sometimes difficult and limited to flow systems which are 
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relativly simple. In actual field conditions, the soils may be 
heterogeneous and anisotropic and boundary conditions very 
complex. These conditions make closed form methods very 
difficult and many simplifications would have to be made. 
However, with the numerical models, these complexities can be 
easily handled and many more check cases can be made in much 
less time. 

Finite Difference Methods 



Many numerical computer models have been developed using finite 
difference methodology. Since the classic work of Peaceman and 
Rachford (1955), the oil industry has sponsored many studies on 
model development. A good historical background of the finite 
difference methods developed in seepage analysis is giving by 
Huyakorn (1983). Worth mentioning are the works of Cooley 
(1971) and Freeze (1971) in the development of other finite 
difference models. These models solve for the free surface and 
treat the satur ated-unsatur ated flow as an immobile air phase. 

Finite Element 

The finite element method for solving subsurface flow has been 
widely used and researched be soils engineers and groundwater 
hydrologists. The program used in this research and discussed 
in more detail in Chapter 3 is a finite element model. 
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Most finite element models are built on the Galerkin method. 
Some earlier finite element models for steady state and 
saturated flow were developed by Desai and Abel (1972), Conner 
and Brebbia (1977) and Zienkiewicz (1977). Finder and Frind 
(1972) used i sopar ametr i c quadrilateral elements with the 
Galerkin" s technique. Neuman (1973, 1974) used a four node 

quadr i 1 ateral element. A good reference for other works in the 
finite element models that have been developed is reported in 
Huyakorn and Finder (1983) and Heijde (1985). 

Some significant works were those of Taylor and Brown (1967) 
where a finite element model was developed for seepage through 
a dam solving for the free surface. The flow was only 
considerd for the saturated zone in this model. Freeze (1971) 
demonstrated that the conditions of the unsaturated zone may 
strongly influence the position of the free surface and results 
obtained not considering the capillary action for this case 
could be in error. Fredlund and Morgenstern (1976) developed 
improved constitutive equations for an unsaturated soil. Lam 
and Fredlund, etal (1984) presented additional theories and 
numerical development for further modeling techniques to 
consider the unsaturated zone in seepage analysis. 
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CHAPTER 



FINITE ELEMENT STUDY 
3.1 Finite Element Program Theory 
Introducti on 

Finite element analysis of seepage problems can be described as 
follows (Huyakon and Finder, 1983): 

(1) A discrete number of i n ter connected nodes to define a 
series of sub-regions or elements. These can take the form of 
triangles, quadr i 1 ateral or rectangular shapes which are so 
spaced to define the geometric boundaries and soil types. 

(2) A matrix expression is developed to relate flow 
variables to each node of the elements. The mathematical 
expression may be of the form of e?ither a variational or 
weighted function. A weighted residual function is used in the 
Galerkin method (Conner and Brebbia, 1977). 

(3) A set of algebraic equations that describe the entire 
global system is applied to assemble the element matrices. 

These equations are derived to consider compatible conditions 
at each node shared be each adjoining element. 

(4) Boundary conditions of flow are incorporated to reflect 
appropriate limits to the global matrix equations. 

(5) The resulting set of simultaneous algebraic equations 
are solved by some alogrithm method. Many algorithyms have 
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been used to e-f f i c i en 1 1 y take advantage of the banded and 
symmetric features of the coefficient matrix (Bathe and Wilson, 

1976) . 

The finite element method has great potential in seepage 
analysis problems, in that the complex geometries of flow 
regions can be easily constructed via nodes and elements. This 
chapter will present the general theories and formulation of 
equations for the model used in this study. 

Program Theory 

A finite element program developed by Dr. Mustusfa Aral of the 
School of Civil Engineering, Georgia Institute of Technology, 
was used to study the flow and variables involved for 
subsurface drains as earlier described. This program was used 
to see effects of different parameters in the analysis performed. 
The program was operated on the Georgia Institute of Technology 
main-frame CDC CYBER computer. A complete user’s guide for the 
program is included as Appendix A-l. Also example input and output 
used in this study is included by Appendix B-l. The formulation 
for the program presented here was summarised from a previous study 
(Aral , Sturm and Fulford, 1981). 

This model only considers the flow condition of the free 
surface being the point where the water level is at atmospheric 
pressure. In other words, this model does not treat the 
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unsaturated zone of the soil. 



For this study the -flow continum 



is modeled to determine the quantity of flow at each drain and 
to determine the free surface location between each of the 
drains. This model also computed the velocities at each 
element for a good check. 



Equations of Flow 

Darcy's Law is used as the basic relationship governing 
groundwater flow. Flow is assumed to be laminar and the 
inertial forces are negligible compared to viscous forces. 
Rewriting Darcy's Law in the cartesian coordinate system: 

v = - K V h (3.1) 



where V is seepage velocity, h is piesometric head and K the 
hydraulic conductivity. In the two dimensional case, equation 
(3.1) can be written as 



3 h 

v = - K ijL , i = 1,2 (3.2) 

9 x 



ii 



where K 



is a conductivity tensor. 



The continuity equation for the flow in an element of soil can 
be expressed 

8 9 

(°v.) = - <n P ) , i = 1,2 (3.3) 

3x 3t 



where T) is porosity of the soil and p is the density of the 
fluid. Assuming the fluid and soil particles are 
incompressible, equation (3.2) and (3.3) can be used for the-? 
following expression in the two-dimensional form: 



a 3 h 

Kii = 0 , i = 1,2 (3.5) 

5x 3 x 



Figure 3.1 can be refered to for notation of the key variables 
used is this study. 
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Boundary Conditions 



For -free surface problems, which is the case for this study, 
the exact location of the phreatic surface is not known at the 
beginning of the solution. Location of the free surface is an 
iterative process using two independent boundary conditions 
available at the free surface. The line of seepage at the free 
surface is the upper streamline of the flow domain. At every 
point along this streamline the pore pressure is at atmospheric 
pressure and therefore the piezometric head equals the elevation 
head. Also the flow across this boundary is zero because of it 
being a streamline. 

Considering the seepage domain as depicted by Figure 3.1, the 
steady state boundary conditions can be described as follows: 

(a) Impervious Boundary (AB) : The program considers this 

boundary as a streamline where velocity to this boundary equals 
zero. 

<b) Constant Head Boundaries (AC) and (EF) : Constant 

hydrostatic pressure is exerted on these boundaries where total 
piezometric head is constant, i.e. equal to the water level in 
a piezometric tube. 

<c) Phreatic Surface (CE) , (FG) , (HI), (JK) and (LM) : This 

boundary is a combination of both cases (a) and (b) above. By 
iteration, described later, both conditions must be satisfied 
independently. 
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(d ) Seepage Faces 
-forces occur at these 
the soil continue, i. 
-flow must be normal t 



(EF) , <GH), <IJ>, <KL) and (MN) 

boundary -faces where water exi 
e. trench drains. Ther-fore str 
o these boundaries. 



: Seepage 

ts through 
eamlines of 
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3.2 Finite Element Fonm.il at i on 
General Equations 

The first step in solving the seepage problem by the finite 
element method is to discretize the flow region into a finite 
number of subregions called elements. Using finite elements, 
equilibrium equations must be made to handle the geometries 
involved. This results in a set of simultaneous algebraic 
equations rather than differential equations. Figure 3.2 
depicts the quadr i 1 ateral elements used for this program and 
shows how four six nodal triangles are formed internal to the 
program. A more detailed description is of finite element 
equations and formulations can be found in Zienkiewicz (1971), 
Desai (1972) and Aral (1974, 1976). 




Figure 3.2 Typical element and nodes 
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A -finite element approximation to equation (3.5) can be written 
through a variational approach. Assuming Dirichlet and Neuman 
boundary conditions, and considering the two dimensional case 
of steady flow, the following equations can be expressed: 



h = h 



(3. 7) 



3 h 



3h 



and V 



n 



K , , a + k 00 a . 

11 ~ nl 22 . ru 

d x o x 



n 



on S 2 (3.7) 



where h is the specified Direchlet boundary condition in terms 
of piezometric head; V n is the velocity in the direction 
normal to the boundary; and & n2 are the direction cosines of 

the normal with respect to Xj and x 2 ; and and S 2 indicate 
respective boundaries of the domain. From equations (3.5) and 
(3.7), the following variational statement can be written: 




+ 




Shdx^ d>: 2 




(3.8) 
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Integrating by parts one obtains 




II 



9h 9 <$h 




8 h 9 £h 




V n 6 hdS 



0 



9 ) 



S 



2 



Since -from earlier discussed boundary conditions as -for this 
problem studied, V n = 0. Thus, equation (3.9) can be -further 
simplified by eliminating the boundary intregrals. 

Assuming equation (3.9) is written -for a single element, matrix 
equations can be generated using an interpolation -function -for 
approx i mati ng h over an element. The interpolation -function 
used in this model is given as 



h 




k 



1 0 ) 



where h^’s are the unknown nodal values o-f the dependent 



variable and N ^ is the interpolating polynomial used to 



approximate h in an element. 
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For a typical six-nodal triangle element these polynomials are 
g i ven as 

N. = £. < 2 £. -i) i = 1,2,3 (3.11) 

111 

and N „ = 4?.?. j =2,3,1 (3.12) 

where ^ > and £3 are area cordinates defined by the 

following relations 








(3. 13) 



where Ai , A 2 and A 3 are the areas of the three subtriangles 
subtended by the point (P) and the corners, the index on (A) 
designating the apposite corner numbers, and A is the area of 
the triangle (123) shown in Figure 3.3. 




Figure 3.3 Area coordinates. 
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To simplify the resultant equations, the origin of the local 
axis which defines the local coordinates of the element, nodes 
are placed at the centroid of the element and the principle 



ani satrapy . 

Substituting equation (3.10) into equation (3.9) the final form 
of the approx i mat i ng equations in an element can be expressed. 
The stationary condition of equation (3.9) is obtained by 
equating the first partial derivatives of the function, with 
respect to the undetermined parameters, hj , to zero: 



axis are inclined angle © in the direction of local 





6 



14) 



This results in a six by six element stiffness matrix, S e , and 



a six by one element load vector F 



e 
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element matrix equations of equation (3.14) is 



The six by six 
wr i tten 



CS e ] £h e > 



v. r j 



<3. 15) 



where 



CEf 3 « 




3 N . 3 N . 

K 22 1 — 3 
3x 2 Sk, 



dXj dx 2 , 



i > j 



(3. 16) 



The same procedure is taken -for all elements. Then the total 
stiffness matrix is arrived at by combining the element 
stiffness matrices using structural assembly techniques 
(Zienkwicz, 1971) 



For a grid of elements, a system of banded simultaneous 
equations results: 



USD Ch > 



jr cr \ 
\ r s 



(3. 17) 



where CB] is the global matrix of cofficients which 
incorporates the properties of the materials and the geometry of 
the elements; Chi- is the vector of the unknown h’s at the 
nodes; and CF> is the load vector. After the introduction of 
appropriate boundary conditions, this system of equations can 
now be solved by computations for the unknowns, h, yielding the 
piezometric head distribution in the solution domain. 
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Free Surface Computation 



The -flow domain must initially be established assuming an 
arbitrary tree surface, i.e. the user initially estimates the 
free surface loacation. All boundary conditions are solved for 
the preset conditions of the problem. The nodes which are not 
fixed are allowed to move, where during the adjustment the 
interior nodes are also moved in order not to have distortion 
of the elements. Then additional iterations are made until all 
boundary conditions and user prescribed error correlations have 
been met for the configuration. 
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3.3 Variables and Dimensional Analysis 

The variables investigated in this study are shown in Figure 
3.1. For the analysis of seepage using the -finite element 
program, variables were changed for each computer run. From 
each run data was generated to compare the changes occuring for 
each variable change. 



Only the two-dimensional flow condition was considered in this 
study. For this flow condition, the building drainage system 
is in one cross-sect i on of the building where each drain 
collects the seepage and collector pipes or drains take the 
drainage by gravity to a sump or some exit drain. This is 
depicted by Figure 3.4. 
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Figure 3.4 Underdrain System For Seepage Removal. 
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The total discharge o-f the system, q T , 
following function: 



can be express as the 



g q. 

where 



= Hr, h , d w , S p , K , Spg, 8 ) (3.18) 

q T = total discharge of all drains, i.e. the sum of 

discharges through each seepage face for the 
model, ft 3 /day/ft 

Hr = depth of saturated soil from the drain 
discharge elevation to the impervious 
boundary, ft 

h = fixed head of above the drains , ft 
d w = drain width, ft 

Sp = drain spacing, ft 

K = hydraulic conductivity of the soil, ft/day 
S FS = difference in elevation of the free surface to 
the water level of each drain, ft 
0 = slope of the impervious boundary with 

horizontal, ft/ft 



Dimensional analysis of equation (3.18) yields: 




(3. 19) 



Since permeability is moved to the left side of equation (3.19) 

these variable is made proportional to the discharge. The 

elevation of the free surface, S FS , is an output variable, so 
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it is not given in equation (3.19. Also since the slope of the 
impervious boundary, 0 , is already dimensionless it can be 
given as a single term in FT/FT. 

The effects of these functional variables and other results are 
given in the form of dimensionless plots in Chapter 4. A 
discussion of the results is given in that chapter. Since 
considerable interest is in quantities of discharge and the 
location of the free surface, this analysis can be used to see 
trends and effects that each variable has on the seepage 
domai n . 
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3.4 Program Ver i f i cat 1 on 

To verify that the computer program used in this study is 
giving valid results, two other methods of solutions were used 
for a selected check case. Another micro computer finite 
element program, SEEP, and an appropriate analytical method 
based on Dupuit’s assumption, was used for the simple 
configuration of a drainage blanket with a single pipe type 
drain (Duncan and Wong, 1985). This is shown in Figure 3.5. 




Figure 3.5 Check case example of a blanket drain with single 
pipe trench drain. 

A grid was generated for the SEEP program to match the 
configuration and boundary conditions as by Figure 3.4. The 
total flow for this check case is calculated by the SEEP 
program as based on the Dupuit assumption of a fixed domain and 
fixed amount of head. Seep also tabulated the flow quantities 

for all nodal points with fixed heads. For this model, the 
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■free surface converged on the -free surface location with an error 
criteria of 0.16 feet. 



The location of the free surface as computed by the Georgia 
Tech model converged to the free surface having an error factor 
of only 0.001 feet, a much more accurate and sensitive 
convergence. However, the two models compared very close as 
far as the location of the free surface. Table 3.1 lists the 
elevations of the free surface nodes for the two models. 
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TABLE 3.1 Elevations of the Free Surface Nodes 



SEEP: VP I 
Model 



Ga Tech 
Model 



x coordinate 



y coordinate 



0 


78. 00 


78.00 


40 


72. 55 


72.73 


70 


67.07 


67 . 57 


80 


64. 27 


65.17 


85 


62. 03 


63.84 


87. 5 


60. 39 


53.73 


90 


58. 00 


58. 00". — drain 

/ 


92 


58.00 


58.00' 


95 


60. 08 


58.32 


97. 5 


61. 16 


62. 15 


100 


61 . 92 


62.52 


106 


63 • OO 


63 . 28 


115 


63.76 


63.93 


130 


64. 26 


64.42 


150 


64.42 


64.62 



1 ocati on 
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To verify that flow quantities were accurate -for the Ga Tech 
model, results were also compared with the results obtained 
-from a Dupuit assumption solution, developed by the Federal 
Highway Admi ni strat i on ’ s, Highway Subdrainage Design Manual, 
(1980). Figure 3.6 was used to calculate the quantity o-f -flow 
to the drain. The -following calculations based on Figure 3.6 
are presented: 



b 1 

H 58’ 

o 



L . 91’ 

l _ 

H 58’ 

o 



1.57 



K ( H - H 0 ) 0.2 (78’ -58’) 



<32 



<32 



2.17 (from Figure 3.5) 



therefore q£ 



0 . 2 ( 20 ) 



2. 17 



1 . 84 CFD 




1.84 



K 



0.2 



9. 22 




(H - H ) 2 

Q 

2 ( L -b) 



( 20 ) 



2(91-1) 



2 . 22 



<3 d 

— = 9.22+2.22 = 11.44 

K 



therfore q d = 0. 2 ( 1 1 . 44) = 2.29 CFD 

The results of this computation is compared with the results 
of the two finite elements in Table 3.2 
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H 

o 

Figure 3.6 Chart for Determining Flow Rate in Symmetrical 
Underdrains (after FHA, 1980) 
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Table 3.2 Discharge to the drain check cases, CFD 





Discharge, CFD 


SEEP: VP I Model 


2. 15 


Ga Tech Model 


2.34 


FHA Nomograph 


2,29 



In summary, the Ga Tech model compares very well with other 
model and with the FHA nomograph -for -flow quantities. The Ga 
Tech model was -found to more error criteria sensitive. The 
•free surface location check also proved to compare very well 
with the SEEP model. 
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CHAPTER 4 



RESULTS AND DISCUSSION 

4.1 Introduction 

This chapter discusses the analysis and results of the two 
dimensional, steady state, free surface finite element study of 
seepage to building underdrain systems. Details presented 
include model set-up such as boudary conditions, grid 
development and basic first assumptions. Output variables 
considered are total discharge and discharge at each drain, 
plus the free surface elevation between the drains. 

In studying sensitivity of the model for the input variables, 
it was found that the distance from the Dirichlet boundary 
face, or constant head boundary, Li, to the first drain was 
critical. Li was found to be critical in terms of output 
charcter i sti cs, i.e. discharge and location of the free surface 
between drains. Therfore, two sets of results are presented for 
two Li distances of 400 feet and 80 feet. 

4.2 Model Set-up 

Figure 4.1 shows the basic grid that was developed for this 
seepage problem. Using 150 elements and 521 nodes, the 
geometry of the mesh was based on a five drain system. The 



CONSTANT HEAD 



UJ 



0-1 

2tr 
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Figure 4.1 Example Mesh Used 



sides and bottoms of the drains were modeled as seepage faces. 
Using the two-dimensional model, the fifth drain is considered 
to be at the building centerline. Assuming symmetry, only 
one-half of the building was analysed. The left most 
boundary, as shown in Figure 4.1, is a Dirichlet boundary 
condition which has the constant head required to allow for the 
seepage to enter the drains. 

Listed below are some basic assumptions for the study made for 
this seepage problem: 

(1) Steady state flow conditions 

(2) Two-dimensional analysis where the right most 
boundary is the building centerline and by symmetry the flows 
are only one-half of the total for the total width of the 

bui lding. 

(3) Free surface nodes move to the computed phreatic 
level, i.e. piezometric head equals atmospheric pressure at the 
upper most nodes 

(4) Partial saturated flow is not considered 

(5) Evaporation is not considered 

(6) Retardation is not considered 

(7) Infiltration is not considered 

(8) The soil is homogeous being either isotropic or 
anisotropic. 
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4 



Variables Used 



For this study, key variables were used in the -finite element 
model to look at their effects on the -flow of seepage for the 
building underdrain system. Two major design parameters 
necessary to model any seepage problem is the total discharge 
of the drain system and the location of the free surface. For 
this study, the following output variables are used: 

q d = discharge or seepage quantity to each drain 
(cubic feet per day per linear foot) 
q T = total discharge of all drains (same units) 

Spg = height of free surface above drain elevation. 

Discharge quantities are in units of flow quantity in time per 
unit length of the total flow continum. Therfore, flow for a 
total system should be multiplied be the length of the 
structure, i.e. in the third dimension. 

Input variables used in this problem to study their effects on 
seepage are listed below: 

Li = distance from the left most constant head 
boundary to the first drain (feet) 

K = hydraulic conductivity (feet /day) 

Kv = hydraulic conductivity in the vertical 
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direction (-feet /day) 

Kh = hydraulic conductivity .in the horizontal 
direction (-Feet /day) 

h = hydraulic head o-F -Flow continum above the 
drain elevation (-Feet) 

Hr = depth of aqui-Fer between the elevation of the 
drain and the impervious boundary (feet) 

Sp = spacing between drains (feet) 

Dw = drain width (feet) 

0 = slope of the impervious boundary ( 7 . ) 

For better reference to the notation of the variables refer to 
Figure 3-1 in Chapter 3. 

For this finite element problem the units of feet and days wer 
used to give results for the output variables in cubic feet, pe 
day and feet per day for velocity. 

Hydraulic conductivity or permeability values used for this 
study were 0.0002, 0.02, 0.2, and 4 feet per day. A soil's 
hydraulic conductivity having a value of 0.0002 feet per day i 
equivalent to a fine grained soil such as a clay. Hydraulic: 
conductivities are tabulated in Table 4.1 to show cross 
reference to soil types and permeabilities in units of 
centimeters per second. 
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Table 4.1 Hydraulic Conductivity Values For Various Soil Types 



Fine 


g r a i n ed 


soi 1 


0 . 0002 


f eet / day 


7x10 


V.m/sec 


Fine 


gr ai ned 


soi 1 


0 . 002 


■feet /day 


7x10 


^cm/sec 


si It 






0 . 02 


f eet /day 


7x10 


^cm/sec 


sand , 


silt 




0 . 2 


f eet /day 


7x10 


-7 

cm /sec 


cl ean 


sand 




4 


f eet/dav 


1 . 4x 10 


-7 

cm/ sec 



The other input variables used were typically realistic values 
to meet typical -field conditions. Listed below are some of 
this values, where for each computor run specific variables 
were changed and the others held constant. The intent was to 
see what effects that each variable change made on the model. 



Li 


= SO 7 , 


400 


7 and 1000 7 


h 


= 10’, 


20 7 


and 40 7 


Hr 


= 20 7 , 


SB 7 


and 100 7 


Sp 


- 15’ , 


25 7 


and 35 7 


Dw 


= 1 7 , 


2 7 and 4 7 


0 


= 07., 


5/C 


107. and 207. 


Kv 


= . 000 


2, 0 


.2, 4 and 1 



Kh = 0.2 

K = 0.0002, 0.02, 4 



59 



4.4 Model Sensitivity 



The model was developed assuming a homogeneous isotropic media, 
except tor a tew cases where anisotropic conditions were studiec. 
The nodes ot the top three rows ot elements were allowed to move 
tor iteration ot the tree surface. In all runs, the maximum 
number ot iterations was six and on an average three to tour were 
found t cd be adequate. The model was made to iterate to an 
accuracy of 0.001 feet during the free surface moves. 

In calibration of free surface models careful attention must be 
made to the moving nodes. Problems exist when the elements become 
less than zero in area or when the elements actual may become so 
distorted that they may become inverted. This was not the case 
•for the mesh used in this study. 

In studying the flow quantities and velocities for the basic 
geometries and variable changes, all results seemed reasonable and 
observed trends appeared physically possible. However large 
variation of output was noticed by the change in distance from the 
constant head boundary to the first drain, Li. The first set of 
results is presented based on the influence length, Li, being 10 
times the maximum head difference, i.e. 10 x 40 7 - 400'. Depth of 

aquifer for this case was held constant at 100 feet. Use of 10 
times the head differential is based on E-iear 7 s (1979) 
recommendat i on . 
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Figure 4.2 depicts the moved -free surface corner nodes for an 
influence length of 400 feet. For comparison, Figure 4.3 shows the? 
moved free surface corner nodes for Li equal to 80 feet and a 
thickness of aquifer, Hr, equal to 58 feet. 
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-20 0 20 40 60 80 100 120 140 160 

HORIZONTAL COORDINATE 

Figure A. 3 Moved Free Surface vs. Original 

Li = 80', Hr = 58', Sp = 25 ' , Dw = 2 ' , K = 0.2 ft/day 



Figure 4.4 and 4.5 are presented to show the dramatic influence 
that the distance from the constant head boundary from the first 
drain, Li, has on total discharge and the maximum free surface 
elevation, respect i vel y . For further reference, the total 
discharge, q , is the sum of the drain discharges and the maximum 
free surface is the highest elevation of water above the drain 
elevation. Lengths between 300 ? and 400'’ appear to result in less 
change in flow quantity and free surface changes. This equates to 
the distance where the constant head boundary away from the drains 
begins to cause the least effect on discharge and free surface 
changes is approx i mat el y 15 to 20 times the head differential. 

This perhaps is good evidence that the actual field condi tons for 
each particular application of groundwater movement should 
carefully be looked at in terms of this distance. 
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Figure 4.4 Distance -from the constant head boundary to the -first 
drain vs. total discharge o-f drains, 
h = 20’, Sp = 15% K = 0.2 -ft/day 
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n o -i FREE SURFACE ELEVATION 
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4.5 Results 



This section will present the computer results in the form of 
dimensionless plots and other nomographs to show the effects 
that each variable change has on the seepage problem. Two 
sets of plots are given where the first set is for the constant 
head boundary distance, Li, equals 400 feet and the depth of 
rock. Hr, equals 100 feet. The second set of plots are 
essentially the same as the first except Li equals 80 feet and 
the Hr equals 58 feet. For variable notation refer to Figure 
3 . 1 . 

Drain Discharge (Li = 400* and Hr = 100") 



Figure 4.6 shows the influence of drain width, Dw, for head 
differentials, h, of 10, 20 and 40 feet on total discharge of 

the drains. This shows that drain width is not particularly 
critical for low head differences. It is reasonable to assume 
that total discharge could approach a constant for drain widths 
greater that four feet and head less than 20 feet. Although 
total discharge increases with drain width, this plot does not 
show what drain width is the most effective or most efficient. 
However, the influence of head can be easily observed 
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q 



T 




Dw 



Figure 4.6 Total drain discharge vs. drain width with head 
di f f erenti al s. 

Li ss 400’ , Hr = 100" , Sp = 25", K - 0.2 ft /day 
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Figure 4.7 depicts the results o-f varying anisotropic soils 
with varying heads. In this analysis, vertical hydraulic 
conductivity varies while horizontal conductivity was held 
constant. Figure 4.7 shows that with increase in vertical 
hydraulic conductivity discharge increases. This appears to 
satisfy theory , that with less resistance flow increases. 




Kv 

Kh 



Figure 4.7 Total drain discharge vs. anisotropy with varying 
head di f f erenti al » 

Li = 400’, Hr = 100’, Sp = 15 ? , Dw = 2’, Kh = 0.2 ft/day 
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Figures 4.8 through 4.12 are graphs which show the discharge 
quantities -for each o-f the -five drains. Figure 4.8 depicts the 
e-f-fects o-f depth of impervious layer, Hr, on the drain 
discharge -for each o-f the drains. This data shows that with 
increased depth discharge increases. It is interesting to note 
that practically no -flow is recorded -for drain 5 when a shallow 
depth o-f impervious 1 ayer exists, i.e. Hr = 20 -feet. 
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Figure 4.3 
Li = 400" , 



Drain discharge vs. thickness o-f aquifer, 
h = 20", Sp = 25", Dw = 2", K = 0.2 -ft/day 
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Figure 4.9 shows the effects of head differential on the 
discharge of each of the drains. It is interesting to see that 
though the first drain is the most critical drain, i . e. a 
large amount of the total flow is collected in this drain, the 
remaining drains collect discharge decreasing linearly from 
the left. This is possiblly due to the influence of the 
constant head boundary distance, Li, being 400 feet away. This 
shows that for the drains, flow quantity increases linearly with 
increase in head. 
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q d DISCHARGE 



F i gur 
Li = 



0 . 18 





2 3 4 5 

DRAINS 



e 4.9 Drain discharge vs. 
400=', Hr = 100’ , Sp = 25’, 



head differential. 

Dw = 2’, K = 0.2 ft /day 
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Figure 4.10 depicts the results of comparing the drain 
discharge versus the spacing between drains. It is of 
interest to see that drain number one for spacing equal to 25 
feet, collects more seepage than the other discharges for the 
other spacings. However, at drains 3 through 5 this changes 
and the discharge for 15 feet spacing is greatest. This will 
be analyzed 1 ater in the dimensional analysis. 
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Figure 4.10 Drain discharge vs. drain 
Li = 400 7 , Hr = 100', h = 20', Dw = 2' 
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spac i ng . 

K = 0.2 ft /day 



.1 



Figure 4.11 shows the effects of anisotropy of the soil 
matrix on drain discharge. It is of interest to see that 
the discharge for Kv/Kh equal one is much greater than the 
other drain discharges for Kv/Kh less than one. 




DRAINS 

Figure 4.11 Drain discharge vs. anisotropy. 

Li = 400’, Hr = 100’, Sp = 15’, Dw = 2, h = 20’ 
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Figure 4.12 depicts data generated -for the condition of a 
sloping impervious layer. The slope of the boundary is given 
as a percentage of Ax/Ay. The depth of the 1 ayer is held 
constant at drain one and the layer slopes clockwise from the 
positive k horizontal axis. 

Though actual quantities of discharge are only changing in very 
small ammounts, the trend is seen that for slopes greater than 
57. the discharge increases in drains 3 through 5. However at 
some slope between 57. and 107 the discharge decreases. It is 
also seen that the discharge of drain five is the least effected 
by any slope change, in terms of discharge quantity. 
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qd DISCHARGE 



0 , 12 




1 2 3 4 5 

DRAINS 




Figure 4.12 Drain discharge vs. slope of the impervious 
1 ayer . 

Li = 400’, Hr = 100’ (first drain), h = 20’, Sp = 25’, Dw 
K = 0.2 ft /day 
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Free Surface Location (Li 



400 ’ and Hr 



1 00 ? ) 



Figure 4.13 shows the influence of the drain width on the 
maximum free surface height, S pc . In all cases this occured 
between drains one and two. This plot shows how the 
relationship is linear and that the height of the free surface 
is not greatly affected by the drain width. 
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D w (ft) 



Figure 4.13 Maximum -free sur-face height vs. drain width with 
varying heads. 

Li = 400% Hr = 100% Sp = 25%, K = 0.2 -ft /day 
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r 2 



Figure 4.14 (depicts the influence 
the maximum -free surface height, 
increased drain spacing causes an 
surface. Therefore one should cons 
possibly deeper drains. 



of drain spacing and head on 
These results show that 
increase in the free 



ider the smaller spacing or 




Figure 4.14 Maximum free surface 
ith varying head, 
i = 400’, Hr = 100’, Dw = 2’, K 



height vs. spacing of drains 
= 0.2 ft/day 
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Figure 4.15 shows the relationship o-f maximum -ft- ee sur-fact; 

height, Sttq , as a -fuction o-f head, h, and anisotropy, Kv/'Kh. 
• max' 

The same trend is seen here as noticed with discharge -for 
anisotropy, that with increased vertical hydraulic 
conductivity the height o-f the -free surface increases. 





Kh 



Figure 4.15 Maximum -free 
varying head. 

Li = 400, Hr = 100', Sp = 



surface height vs. 
25 ? , Dw = 2 7 



anisotropy with 
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In Figures 4.16 through 4. 20, the -Free surface corner nodes are 
plotted -for the diFFerent cases o-F analysis. To see the 
effects o-F the Free surface changes versus the variable changes 
only nodes between the drains are plotted. In the sence oF 
design, the maximum height oF the Free surFace would be oF the 
most concern, which is in between the First two drains. 

However, it is oF interest to see the total Free surFace From a 
standpoint oF total cross section oF the building. 

In Figure 4.16 the eFFect oF thickness oF aquiFer on location 
oF the Free surFace is shown. It is oF interest to see, as in 
the case For discharge shown by Figure 4.8, that For shallow depths 
oF impervious layer, Hr = 20 7 , that very little rise in the 
Free surFace occurs in between drains Four and Five. 

Figure 4.17 shows the eFFects oF head diFFerential on the 
location oF the Free surFace nodes. It is oF interest to see 
For cases oF head equal to 20 and 10 Feet the Free surFace 
nodes are approaching gentle mounds oF water in between drains. 
However For the head equal to 40 Feet the Free surFace peaks 
rather abruptly between the First two drains but resumes the 
gentle Fall in Free surFace between the other drains. 

Figure 4.18 shows the eFFects oF drain spacing on the location 
oF the Free surFace nodes. This shows as spacing is increased 
the Free surFace increases. 
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For Figure 4.19, the effects of the sloping aquifer on location 
of the free surface nodes is shown. The same trend is seen 
here as in discharge, shown by Figure 4.12, that for slopes 
greater than 10/C the free surface location drops. The change 
in free surface nodes is very small for the cases of slope 
equal to 57. and 107.. 
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Figure 4.16 Free Surface corner nodes at drains with effects 
impervious boundary layer. 

Li = 400' , h = 20' , Dw = 2 \„ K = 0.2 ft/day 
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Dimensionsal Analysis (Li 



400 ? and Hr 



100 ? ) 



Figures 4.21 through 4.24 are plots of the data generated by a 
dimensional analysis as discussed in Chapter 3. The discharge 
values given are unit discharge quantities and each case is -for 
the particular variables noted in each -figure. Unit discharge 
is defined as the sum o-f the total drain discharges -for each 
o-f the drains per unit width in the z axis or the third 
di mensi on . 



In Figure 4.21 the effect of drain width is analyzed. This 
plot shows that with increased drain width discharge increases, 
but with increased head the discharge may increase at a much 
faster rate. From this analysis it would appear that the 
greatest increases of flow rate would be for the drain widths 
between one and two feet wide. The following analysis is 
per f ormed : 



Drain Widths 



head 1 



10 0.034 0.053 0.065 0.073 

increase - 0 .019 0 . 0 12 0 . 009 

20 0.035 0.055 0.070 0.086 

increase *- 0.02 0.015 0.016 

40 0.036 0.057 0.075 0.09 



increase •- 0.021 



0 .018 

88 



0 .015 



This analysis shows the increases o-f the -flow -function -for the 
di-f-ferent head increases as shown. The maximum increase is 
seen -for the values between drain widths o-f one and two. 
Ther-fore the most effective increase in -flow occurs -for drain 
widths between 1 and 2 -feet. Another relationship could be 
stated, that for Dw/h ratios less than 0.2 the flow quantity 
increases at a rapid rate. 
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In Figure 4.22 the effects of thickness of aquifer is shown for 
the dimensional analysis. This plot concludes that for Hr/ h 
ratios greater than 2 or 3 the net increase in flow is 
approaching a linear rel at i onshi p . This states that with 
increasing thickness of the aquifer the dicharge is expected to 
increase very near a linear increase and for very shallow 
thicknesses the flow becomes very small. 




Figure 4.22 Dimensional Analysis, Depths of impervious layer. 
Li = 400’ , Sp = 25’, Dw = 2’, K = 0.2 ft/day 
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The 



The effects o-f drain spacing is presented in Figure 4.23. 
flow function for this analysis is treated differently for the 
spacing changes. As the spacing was increased for each of the 
computer runs the flow domain increased accordingly to 
accommodate the five drains at the appropriate spacing. 
Therfore to gain understanding of the effects of flow 
relationships with respect to drain spacing, each term was 
normalized by the corr espondi ng spacing between the drains for 
each head increase. This analysis shows for decreased spacing 
discharge increases and as head increases so the discharge. 
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Hr /h 



Figure 4.23 Dimensional Analysis, Drain spacing. 

Li = 400, Hr = 100’, Dw = 2’, K = 0.2 ft/day 

Figure 4.24 shows the dimensional analysis -for ef-fects of the 
sloping aquifer. For this particular flow matrix as modeled, 
a definate break in flow increase is noticed. The peak 
discharge function is at slope, 0 equal 0.05. 
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Figure 4.24 Dimensional Analysis, Slope o-f aqui-fer. 

Li = 400’ , Hr = 100’ (-first drain) , Sp = 25’, Dw = 2’, K = 
f t/day 
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Results For Li 



SO 7 and Hr 



58’ 



The -following plots and nomographs are the results for the 
computor runs where distance from the constant head boundary to 
the first drain, Li, equals 80 feet and the thickness of the 
aquifer equals 58 feet. The plots are essentially the same as 
for the previous case except these variables are held 
constant . 

Drain Discharge ( Li = 80’ and Hr - 58’ ) 



Figure 4.25 shows the effects of drain width on total drain 
discharge. The results are generally the same as shown in 
Figure 4.6 except discharge quantities range from 3 to 4 times 
greater . 
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+ Head 
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Figure 4.25 Total drain discharge vs. drain width with 
head . 

Li = 80 ? , Hr = 58' , Sp = 25% K = 0.2 -ft/day 



varyi ng 
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The results of anisotropy of the soil continum is shown by 
Figure 4.26 for the short Li distance of 80 feet. The trend is 
slightly different than for the case for Li equal 400 feet as 
shown by Figure 4.7. For the largest head value of 40 feet the 
discharge decreases just slightly with decrease in vertical 
hydraulic conductivity, but for a Kv/Kh ratio less than 0.1 the 
discharge decrease is very dramatic. However for the low head 
value of 10 feet the general decreasing trend is prevelant as 
the case for Figure 4.7. 



For this case the flow quantities have increased 
4 to 5 times. This shows the large effects that 
of the constant head boundary has on the seepage 



approx i matel y 
the distance 
quant i t i es. 
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Figure 4.26 Total Drain Discharge vs. Anisotropy with varying heads 
Li = 80', Hr = 58', Sp = 15', Dw = 2', Kh = 0.2 ft/day 



Figures 4.27 and 4.28 depict the same general trend of maximum 
flow at drain one and a general decrease in flow for each 
consectutive drain is noticed for effects of thickness of 
aquifier and head respectively. Flow quantities are again 3 to 
4 times greater for this case than for the earlier case as 
shown by Figures 4.8 and 4.9. For this case of shorter 
distance to the constant head boundary equal 80 feet, the 
mid-peak discharges at the third drain is not noticed as in 
Figure 4.9. In Figure 4.28 it is of interest to see that the 
largest head value equal 40 feet has such a large effect on the 
first drain collecting a large portion of the total discharge 
of the system. 
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q d DISCHARGE 



0.5 




DRAINS 



Figure 4 
Li = 80 ? 



27 Drain discharge vs. 
h = 20 J , Sp = 25’ , K = 



depth of aquifer. 
0.2 ft/day, Dw = 2 ’ 
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qd DISCHARGE 



0 .? 




DRAINS 

Figure 4.28 Drain discharge vs. head di -f -f erent i al . 

Li = 80 7 , Hr = 58 7 , Sp = 25 7 , Dw = 2 7 , K = 0.2 -ft/day 
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The effects of drain spacing for distance of the constant head 
boundary to the first drain, Li, equal 80 feet, is shown by 
Figure 4.29. The same trends are seen for this case as that 
shown by Figure 4.10. Flow increase is somewhat less in this 
case compared to the earlier case of Li equal 400 feet. Flow 
quantities appear to range from 2 to 3 times as great. 
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SPACING 




15' 


1 


SPACING 




25’ 


1 


SPACING 




35' 



2 3 4 

DRAINS 



Figure 4.29 Drain discharge vs, 
Li = 80’, Hr =58’, h = 20’, Dw 



drain spacing. 

= 2’, K = 0.2 ft/day 
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In Figure 4.30 the effects of slope of the impervious boundary, 
® , is shown for drain discharge. A significant change is 
noticed for this plot when compared to the earlier plot in 
Figure 4.12 for sloping aquifer. For Li equal to 400 feet, the 
peak discharge occured for a slope between 5 and 10 percent. 

In this case, where Li equals 80 feet, the peak flow is for the 
case of the steepest slope equal to 20 percent. For drain one, 
the flow is 8 times greater than the flow for drain one as 
shown by Figure 4.12, Li = 400 feet. The flow is 5 times 
greater in drain five in this case than the flow in drain five 
for the case of Li = 400 feet. 
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qd DISCHARGE 



0.6 




DRAINS 



Figure 4.30 Drain discharge vs. slope of the impervious 
boundary . 

Li = 80S Hr = 58 7 (first drain), h = 20 7 , Sp = 25 7 , Dw = 2 7 , 
K = 0.2 ft/day 
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Free Sur-face Location < Li 



SO’ and Hr 



58’ ) 



Figures 4.31 through 4.36 show the results o-f the -free surface 
location plots for the case of Li equal 80 feet and Hr equal 58 
feet. The same general trends are shown for this plots as 
shown by the earlier case of Li equal 400 feet and Hr equal 100 
feet. In general , the free surface has moved up approximately 
three times in height. An exception to the earlier case of is 
shown in Figure 4.36. This figure shows the location of the 
free surface nodes with the effects of slope of the impervious 
boundary. It is somewhat different in this case than in the 
earlier case as shown in Figure 4.19 for Li equal 400 feet. In 
this case the free surface continues to increase in height as 
the slope of the impervious boundary increases. 
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Figure 4 
Li = 80 ? 



31 Maximum -free surface height vs. drain 
Hr = 58\, Sp = 25’ , K = 0.2 -ft/day 



Head 

10 ' 

+ Head 
20 ' 

* Head 

40 ' 



width. 
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HEAD (FT) 



Figure 4.32 Maximum free surface height v^=>. frain spacing 
with head varying. 

Li = 80 ? , Hr = 58 7 , Dw = 2 ? , K = 0.2 ft/day 
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Depth of impervious boundary. 
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Dimensional Analysis ( Li = 80'' and Hr = 58’ ) 

The effects of drain width as shown by the dimensional analysis 
is shown by Figure 4.37. The general trend of the discharge 
increasing with drain width is seen in the data. When compared 
to the earlier case for Li equal to 400 feet, it is seen that 
total discharge is greatly effected, by a factor of three 
ti mes. H , 4 t ? 





Figure 4.37 Dimensional Analysis 
Li = 80’, Hr = 58’, Bp = 25’, K = 



drain width 
0.2 ft /day 
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Figure 4.38 shows the dimensional analysis results -for 
thichness of the aquifer. In this case, Li equal 80 -feet, the 
relationship is more linear than in the case o-f Li equal 400 
•feet as shown in Figure 4.22. This results in the conclusion 
that the effects o-f Hr has a pronounce influence on the -flow 
quantities when coupled with the distance Li. For a shorter 
distance o-f Li, the effects o-f depths o-f Hr is more crictical 
and more of a linear relationship. However when the distance 
from the constant head boundary is larger, as in the earlier- 
case, the effects of Hr are not as effective when Hr is very 
smal 1 . 
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Hr 

h 



Figure 4 
Li = 80’ 



38 Dimensional Analysis, depth o-f aqui-fer. 
Sp = 25’, Dw = 2', K = 0.2 -ft/day 
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Figure 4.39 shows the results of the dimensional analysis for 
drain spacing. Again the ordinate is normalised with respect 
to the area subtended by each flow continum for each case of 
spacing. Flow quantity would equate to an increase of 
approximately 5 to 6 times from this case of Li equal 80 feet 
versus Li equal 400 as shown in Figure 4.23. 







Figure 4.39 Dimensional Analysis, drain spacing. 
Li = 80’, Hr = 58’, Dw = 2’, K = 0.2 ft/day 
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Figure 4.40 shows the dimensional analysis results tor the 
sloping impervious boundary. For this case, Li equal 80 teet, 
as depicted in the earlier plots tor sloping aquiter tor drain 
dischrge and tree surtace location, the trend is seen that tor 
and increase in slope the discharge increases. 




© 



Figure 4.40 Dimesional Analysis, sloping impervious boundary. 
Li = 80’, Hr = 58’, Sp = 25’, Dw = 2, K = 0.2 tt/day 
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4.6 Summary 



This chapter has included the results of many computer runs -for 
the analysis of a building underdrain system. The importance 
of model sensitivity must be stressed, in that key variables may 
effect the results very significantly. In this model, it was 
found that the most crictical variables affecting the total 
discharges of the drains and location of the free surface was 
the distance from the constant head boundary to the drain 
system, Li, and thickness of the aquifer, Hr. 

The results obtained appear to be reasonable and within limits 
of possible field occurances. However, using these results for 
quantifying flows and the free surface locations requires a 
very careful study and background investigation of all 
geometric configurations and soil conditions. It is intended 
that this study has presented the results in a manner that one 
can observe the effects of the variables studied on such a 
underdrain system. 

In summary, the following observations are discussed for the 
variables used in this analysis. 

Constant Head Boundary Distance, Li 

As distance from the constant head boundary to the system 
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increases the flow quantities increase for the drain system. 

In terms of a correlation with hydraulic head differential 
above the drain, as used by other authors, this model study 
shows that for a distance of approx i matel y 15 times the head, 
discharge appears to decrease. In terms of gradient, this 
equates to a 1 over 15 ratio or 0.067 gradient. A slightly 
higher distance was observed in the analysis of this influence 
to location of the free surface. 

Depth of Aquifer, Hr 

The results show that with increased aquifer thickness the 
drain discharge increases, increasing very near to a linear 
relationship. It was also observed that for the very shallow 
thicknesses, the last drains collected very little flow. The 
height of the free surface is affected in the same fashion. 

With increased thickness of aquifer the height of the free 
surface increases between drains. 

Drain Spacing, Sp 

Drain spacing effects on discharge was greatly affected by the 

distance of the system from the constant head boundary. A 

substantial increase of 5 to 6 times the discharge was observed 

for the two distances used in the results presented. It was 

observed that as head increases, discharge increases for all 

drain spacings; and as drain spacing increases, discharge 
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decreases. It was also observed that -from the dimensional 
analysis, the most efficient drain spacing is for a spacing 
over head ratio of less than 2. The location of the free surface 
was observed to increase between drains as drain spacing is 
1 ncr eased . 

Drain Width, Dw 

With increased drain width the discharge increases for all 
other variables constant. However, there must be a trade off 
for efficiency of the complete system. It appears that a drain 
width between 1 to 2 feet is the most effective width. This is 
also a plus from a standpoint that lateral seepage can be 
easily drained into thinner trenches or deeper constructed 
dr ai ns. 

Slope of Impervious Boundary, 0 

For the distance of the constant head boundary to the drain 
system equal to 80 feet, the drain discharge was observed to 
increase with increased slope. However, for the distance, Li 
equal to 400 feet, the maximum discharge occurs at a slope 
between 5 to 10 percent. This occurance is not fully 
understood and perhaps should be further studied to develope 
some reasoning. This requires that careful consi der at i on of 
specific boundary geometries be used in model studies for seepage 



anal ysi s. 
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CHAPTER 5 



CONCLUSS I ON 



In this day and age, more and more sites are chosen -for 
buildings and engineered structures that are less desirable to 
work with -from a standpoint o-f site conditions, poor soils and 
water problems. A less desirable site may easily be where the 
groundwater table is very near to the surface. Therfore the 
building’s design must incorporate some thought into how the 
groundwater will be handled and what effects it may have on the 
structure. 

Engineers faced with the analysis of seepage for buildings to 
be built below the groundwater table must consider such 
conditions of uplift pressures, critical velocities, seepage 
control of the water itself or particle migration. Hopefully 
the design will enhance stability of the structure and prolong 
the usefulness and life of the structure. 

In performing the literature review on this topic, it was found 

that many authors and researchers are more and more relying on 

and developing computor models to perform the analysis of 

complex seepage problems. With the use of general seepage 

theory and mathematical representat i on of the seepage theories, 
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the computer programs can quickly and ef f i ci entl y aid the 
endgineer in studying the flow conditions. As in this study, 
the analysis of one particular seepage problem using a finite 
element program made it possible to study many key variables 
and develope many rel at i onshi ps . 

For this study, underdrain systems for buildings built below 
the groundwater table were researched. By using a finite 
element program, and by developing a mesh for this particular 
flow matrix, many variables of the analysis were studied. It 
was found that the model could easily be used to efficiently 
and accurately study these variables. Of major concern, the 
di charge of the system and location of the free surface was 
studied in detail versus the various varible analysis. 

This research has attempted to present general seepage theory 
and useful criteria necessary to analyzing this particular 
seepage problem. It has also attempted to use state of the art 
techniques, i.e. finite element models, to analyze the 
variables of the underdrain system design. In conclusion, the 
author must urge that in using such models that the basic 
assumptions of the model be thoroughly understood due to the 
somewhat si mpl i f i cat i ons of the some model operations, and that 
the engineer obtain as much information and data for the soil 
parameters and fluid boundary conditions which may exist. 



120 



REFERENCES 



Aral , M. M. , "Finite Element Solutions of Selected Partial 
Di -f-ferent i al Equations: FEMAC Computor Program, " Publ. No. 

28, Middle East Technical University, 1974. 

Aral, M„ M. , and Isilgan, "Seepage Through Earth Dams: A 

Finite Element Analysis," METU Journal, 1973. 

Aral, M. M. and Masilia, M. L. , "Unsteady Seepage Analysis of 
Wallace Dam," ASCE, Hydraulic Engineering Journal, Vol . 109, 

1983. 

Aral, M. M. ; Sturm, T. W. ; and Fulford, J. M. , "Analysis of the 
Development of Shallow Groundwater Supplies by Pumping From 
Ponds," School of Civil Engineering: Environmental Resources 

Center, Georgia Institute of Technology, ERC-02-81, 1981. 

Bathe, K. J. and Khoshgottaar , M. R. , "Finite Element Free 
Surface Seepage Analysis Without Mesh Iteration, " Inter. 
Journal for Numerical and Analytical Methods in Geomechan i cs, 
Vol . 3, 1979 

Bear, Jacob, Hydraulics of Groundwater , Mcraw— Hill 
International Book Company, New York, 1979. 

Bertram, G. E. , "An Experimental Investigation of Protective 
Filters," Publications of the Graduate School of Engineering, 
Harvard University, No. 267, 1940. 

Bouwer, H. , Groundwater Hydrology , Mcraw-Hi 1 1 Book Co., Inc., 
New York, N.Y., 1978. 



Brooks, R. H. , and Corey, A. T. , "Hydraulic Properties of 
Porous Media," Colordao State University, Hydrology Papers 
No. 3, Fort Collins, Colorado, 1964. 

Carman, P. C. , Flow of Gases Through Porous Media , Academic, 
New York, N. Y. , 1956. 

Casagrande, A., "Control of Seepage Through Foundations and 
Abutments of Dams," Geotechnique, Vol. 11, Sept. 1961. 

Casagrande, A., "Seepage Through Dams," Harvard University 
Publication 209, reprinted from Journal of the New England 
Water Works Association, June 1937. 



122 



Casagrande, A. "The Role of the Calculated Risk in Earthwork 
and Foundation Engineering," the Second Terazaghi Lecture 
presented at the A.S.C.E. Annual Meeting in New York, October 
21, 1964, Journal o-f the Soil Mechanics and Foundations 

Division, Proceedings, A.S.C.E., Vol. 91, No. SM4, July 1965, 
Part 1. 

Cedergreen, H. R. , Seepage, Drainage and Flow Nets , John Wiley 
and Sons, Inc., New York, N.Y., 1977 

Childs, F. C. , and Col 1 i s-George, N. , "The Permeability of 
Porous Media," Royal Society of London, Proceedings, Series 
A, Vol. 201, 1950. 

Connor, J. J. and Brebbia, C. A., Finite Element Techniques For 
Fluid Flow , But ter worth & Co. , 1977. 

Cooley, R. L. , "A Finite Difference Method For Unsteady Flow in 
Variable Saturated Porous Media: Application to a Single 
Pumping Well," Water Resource Research, 7(6), 1971. 

Desai , C. S. and Apel , J. A., Introduction to the Finite 
Element Method , Van Nostrand Reinholdt Company, New York, 
N.Y., 1972. 

Duncan, J. M. and Wong, K. S. , "A Computer Program For Seepage 
Analysis of Saturated Free Surface or Confined Flow," 

Virginia Polytechnic Institute, Dept of Civil Engineering, 
Geotechnical, 2nd Edition, 1985. 

Federal Highway Administration, Highway Subdrainaqe Design , 

U.S. Dept, of Transporation, Report No. FHWA-TS-80-224 , 

1980. 

Ferguson, H. , and Gardner, W. R. , "Diffusion Theory Applied to 
Water Flow Data Obtained Using Gamma Ray Absorption," Soil 
Science Society of America, Proceedings, Vol. 27, 1970. 

Freeze, R. A., "Influence of the Unsaturated Flow Domain 

Seepage Through Earth Dams," Water Resources Research, Vol. 

7, No. 2, 1971. 

Freeze, R. A., "Three-di mensi onal , Transient, 

Satur ated-Unsatur ated Flow in a Groundwater Basin, Water 
Resource Research, 7, 1971. 



123 



Freeze, R. A. and Witherspoon, P. A., "Theoretical Analysis of 
Regional Groundwater Flow:l. Analytical and Numerical 
Solutions to the Mathematical Model." Water Resources 
Research, 2, 1966. 

Gardner, W. R. , "Some Steady State Solutions of the Unsaturated 
Moisture Flow Equation With Application to Evaporation From A 
Water Table," Soil Science Society of America, Vol . 85, 

1958. 



Harr, M. E. , Groundwater and Seepage, McGraw-Hill Book Co., 
Inc., New York, N.Y., 1962. 

Heijde, P. etal , Groundwater Management: The Use of Numerical 

Model s , American Geophysical Union Water Resources Monograph 
5, 1985. 

Huyakorn, P. S. , and Pinder, G. F. , "Computional Methods in 
Subsurface Flow," Academic Press, New York, N.Y., 1983. 

Irmay, S. , "On the Hydraulic Conductivity of Unsaturated 
Soils," Transactions of AGU, No. 1, 1954. 

Kay, B. D. and Law, P. E. , "Measurement of the Total Suction of 
Soil by a Thermistor Pssssssychr ometer , " Soil Science Society 
of America, Proceedings, Vol. 34, 1970. 

Kozeny, J., "Ueber Kapillare Leitung des Wassers im Boden, " 
Wien, Akad. Wiss., Vol. 136, Part 2a, 1927. 

Lam, L. and Fredlund, D. G. , "Saturated— Unsaturated Transient 
Finite Element Seepage Model For Geotechnical Engineering," 
Depart.. Civil Eng., Univ. of Saskatchewan , Saskatoon, 

Canada, S7N0W0, Finite Elements in Water Resources, 
Proceedings of the 5th Inter. Confrence, Burlington, 

Vermont , 1 984 . 

Lambe, T. W. and Whitman, R. V., Soil Mechanics , John Wiley and 
Sons, 1969. 

Leonards, G. A., Foundation Engineering , McGraw-Hill Book Co., 
New York, N.Y., 1962. 

Masha, M. L. , "Numerical Modeling of Satur ated-Unsatur ated 
Fluid Flow Through Porous Media," thesis presented to GA 
Tech, Atlanta, Ga. , 1980. 



124 



Mavis, F. T. and T. F*. Tsui, "Percolation and Capillary 

Movement o-f Water Through Sand Prisms," Bulletin 18, Univ. of 
Iowa, Studies in Engineering, Iowa City, 19939. 

Meinzer , 0. E. , Hydrol ogy , Dover, New York, 1942. 

Milligan, V., "Field Measurement of Permeability in Soil and 
Rock," State of the Art Paper, Proceedings, ASCE Conference 
on in Situ Measurement of Soil Properties, Vol . II, N. . 

State Univ., 1975. 

Musk at, M. , Flow of Homogeneons Fluids Through Porous Med ia, 
McGraw-Hi 1 1 Book Co. , New York , N . Y . , 1937. 

Moulton, L. K. , "Design of Subsurface Drainage Systems for the 
Control of Groundwater, " Presented at the 58th Annual Meeting 
of the Transportati on Research Board, Washington, D.C., 

1979. 

Department of the Navy, Naval Facilities Engineering Command, 
Soil Mechanics , Design Manual 7.1, Alexandria, Virginia, 

1982. 

F'eaceman, D. W. and Rachford, H. H. , Jr., "The Numerical 

Solution of Paraholic and Elliptic Differential Equations," 
Journal of Society of Industrial Applied Mathematics, 3, 

1955. 



F'inder, G. F. and Gray, W. G. , Finite Element Simulation in 
Surface and Subsurface Hydrology, Accademic Press, New York, 
N.Y., 1977. 



F'inder, G. F. and Frind, E. 0. , "Applications of Galerkin’s 
Procedure to Aquifer Analysis," Water Resources Research, 

Vol . 8, No. 1, 1972. 

F'dubar inova-Kochi na, F’. Ya. , Theory of Groundwater Movement , 
Gostekhi zdat , Moscow, 1952; English translation by R.J.M. 
DeWiest, Princeton Univ. Press, Princeton, N.J., 1962. 

Reeves, M. and Dungund, J. 0., "Water Movement Through 

Saturated-Unsaturated Porous Media: A Finite Element Gal er kin 
Model," Oak Ridge National Laboratory Report, No. 4927, 

1973. 



125 



Ran ki lor, P. R. , Membranes in Ground Engineering, Man stock 
Geotechnical Consultancy Services, Ltd . , John - Wiley & Sons, 
New York, N.Y., 1981. 

Sherard, J. L. , Decker, R. S. , and Ryker, N. L. , "Piping in 
Earth Dams of Dispersive Clay," Proceedings, A.S.C.E. 
Specialty Conference on the Performance of Earth and 
Ear th-Suppor ted Structures, Purdue University, June, 1972. 

Si 1 i n— Bekchur i n , A. I.., Dynamics of Ground Water , Moscow 
Izdat., Moscow Univ., 1958. 

Sowers, G. F., Introductory Soil Mechanics and Foundations , 4th 
Ed. MacMi Ilian Publishing Co., Inc., New York, 1979. 

Taylor, D. W. , Fundamentals of Soil Mechanics , Wiley, New York, 
N.Y., 199948. 

Terzaghi, K. , Theoretical Soil Mechanics , John Wiley ?< Sons, 

New York, 1943. 

Terzaghi, K. and Peck, R. B. , Soil Mechanics in Engineering 
Practice , Wiley, New York, N.Y., 1948. 

Than i kachal am , V. R. , and Sakthi vadi vel , R. , "Rational Design 
Criteria for Protective Filters," Canadian Geotechnical 
Journal , Ottawa, Vol 11, No. 2, 1974. 

U.S. Army Corps of Engineers, "Drainage and Erosion 

Control —Subsurf ace Drainage Facilities for Airfields," Part 
XIII, Chapter 2, Engineering Manual , Military Construct i on , 
Washington, D.C., 1955. 

U.S. Army Corps of Engineers, Waterways Experiment Station, 
"investigation of Filter Requirements for Underdrains, " 
Technical Memorandum No. 183-1, 1941. 

U.S. Bureau of Reclamation, Design of Small Dams , U.S. 

Government Printing Office, Washington, D.C., 2nd Ed. 1973. 

Zienkiewicz, 0. C. , The Finite Element Method in Engineering 
Science, McGraw-Hill Publishing Co., 1971. 

Zienkiewicz, 0. C. , The Finite Element Method , Third Edition, 
McGraw-Hill Book Company (UK) Limited, London, 1977. 



126 



APPENDIX A— 1 



USER'S MANUAL FOR 
ARAL SEEPAGE PROGRAM 



I. INTRODUCTION: 

This seepage program is a -finite element model f or solving 
unconfined or confined two dimensional or axisymetric seepage 
problems. For unconfined flow analysis the model iterates by 
moving the free surface nodes to the phreatic level. 

In general, the program allows f<br the porous soil mass to be 
divided into discrete elements and nodes using either a vertical 
plane model (two-dimensional ,x-y) or axisymmetric grid structure. 
Then applying known boundary conditions, the program 
simultaneously solves for each element, the continuity equation 
which satisfies Darcy's Law. 

Figure 1 shows a very simple case of a two dimensional unconfined 
steady state condition which can serve as an example. In this 
example the free surface location is unknown and the user must 
begin the process by i ni t i al 1 i 2 i ng the input by approx i mat i ng the 
location. Then the program will iterate to the free surface level 
until the free surface converges to satisfy the equations of 
f 1 ow . 
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Figure 1 



Example Grid -for 



Uncon-f i ned 



Flow Conti num. 



> 
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II. PROGRAM EXECUTION AND ANALYSIS: 

The user has the option of performing a two dimensional problem or 
an ax i symmetr i c problem. The two dimensional case simply shown by 
figure 1 will serve as an example to explain the parameters. An 
axisymetric application might be a well source problem looking at 
stratified soils having varying permeabilities at different 
elevations through the porous media. 

III. BOUNDARY CONDITIONS OF FLOW: 

The potential, of flow must be defined as some directional 



A Dirichlet boundary condition is where the piezometric head, h in- 
constant, or 



is solved for at a particular point. This boundary condition 
occurs whenever the flow domain is adjacent to an open body of 
fluid. Since the piezometric head is constant at the open body of 



function 




h 




water , the piezometric head of the porous media at the interface 
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will be constant. This boundary condition denotes a non- zero 
•flux or -flow is to be considered possible. 

Another boundary condition -for -flow potential is the Neuman 
boundaries. This is the case where -flow cannot move through the 
porous media at a specific point defined, i.e. an impervious lay or 
of stratum. 

IV. FINITE ELEMENT GRID NOTATION: 



The porous soil mass to be analysed must 
grid of elements and nodes. The program 
i sopar ametr i c element, shown by figure 2 
element containing four corner nodes and 



be divided into a 

uses any shape second-order 

as a quadr al ater al 

four midpoint nodes. 




Figure 2 Element and Node Notation. 



The nodes must be numbered -from some beginning corner, preferably 
the number one node should be at the bottom left corner of the 
mesh. Then node numbering increases to the end of each column, 
beginning at the bottom of each column throughout the mesh. One 
can see the numbering scheme by the examples presented. Elements 
are numbered the same as the nodes. 

The program automat i cal 1 y generates nodes at equal spacing between 
corner nodes which are defined by the user. This feature 
requires that careful attention be made to the node numbering scheme 
once corner nodes have been established. This feature is good for 
creating finer mesh grids in particular areas of interest. 

Base elements are used to create stationary reference elements 
typically at the bottom of the mesh or away from the free surface 
elements. These elements are also used to define changes in 
hydraulic conductivity in the x and y di recti ons (Kx and Ky) . 

The free surface or phreatic surface of a media is defined by the 
user as a first appr ox i mat i on of it's actual location. In the case 
of a cutoff wall, shown in Figure 3, there is no free surface and 
the seepage model solves for the flow char acter i st i cs of each 
element and node. However, in the case of a dam as shown in 
Figure 3, the free surface is unknown and the program will solve 
for the phreatic surface where the fluid pressure is at 



atmospheric pressure. In con-fined aquifers the water pressures 
would possible be higher than atmospheric pressure. 

Seepage -faces are de-fined by the user to denote a boundary or 
where seepage is. to be allowed. Here flow is qun^ati-fied and th 
model solves -for the velocity and discharge at this -face. 




i 



CUTOFF WALL 



DAM CROSS-SECTION 



Figure 3 Example o-f con-fined -flow and unconfined -flow. 



Permeability o-f the porous media o-ften depends on the direction 
of flow. For the condition of soil matrix being anisotropic the 
velocities and gradients are generally not parallel. For this 
condition the velocities can be written in the Darcy Law form: 



Vx 



-Kx x 



^ h 
<3 x 



Vy = 



Vz 



-Kyy 



— . I*:’’ 

I *. A- A- 



<^h 

dy 

<^h 



Normally the hydraulic conductivity of teh media is similar in the 
y direction. Therfore treatment of this normal condition requires 
use of only two direcetions. For this simplification, the program 
uses the notation of the two dimensional case. Figure 4 shows the 
notation that should be used for describing the hydraulic 
conductivity conditions. 



y 




Figure 4 Hydraulic Conductivity Notation. 
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INPUT VARIABLES AND DATA INPUT 



Input can be given in unformated notation where each data i 
seperated by a comma or space. Careful attention, however, 
be made to use the apppropriate number input data per line 
described by this guide. 



nput is 
must 
as 



Title and Type Problem Card. 

TITLE The title can be up to 72 characters. 

DIMEN Flow problem e problem is a two dimensional analysis 
input 2 and for an axisymetric problem input 1. 

Plane Flow Problem 2 

Axisymetric Problem 1 

FERR Having performed a sensitivity of the model or decided 
on an error critiria, input the error required. 



Control cards 

NNPC Total number of corner nodal points. 

These nodes are used to define corners of elements 
which are defined by input coordinates. Nodes that 
are not listed will be generated between two corner 
nodes for each vertical column or row deliniation. 

NELEMC Total number of base elements. 

Base elements are typically at the bottom of a series 
of elements that change the hydraulic conductivity of 
the elements above these base elements. These elements 
are established as reference elements for the program 
to describe permeability of each element. 



NPBDC Total number of Neuman Boundary Faces. 

A Neuman face is defined as a line of nodes that may 
create a barrior or impervious wall for the flow. If 
two corner nodes define a line of several nodes that 
makes an impervious barrier, this is one Neuman 
Boundary face. 

NPFS Total number of corner nodes on the free surface. 

The user must specify particular corner nodes which 
reflect the first appr ox i mat i on of the unknown free 
surface. The program uses these nodes and the midpoint 
nodes to converge to the free surface. 



ITGIV Number o-f iterations desired. 

This is a safety feature for possible input errors and 
computor time limitation. 

NPBDC Total number of Dirichlet Boundary faces. 

Flow is allowed at these faces. Do not count the mid 
point nodes as part of the face. 

NSF'FAC Total number of seepage faces with Dirichlet 

boundaries. For each face defined by two corner no;es 
for the NPBDC input variable a NSPFAC value should be 
given. As in 

NPBDC = 2 NSPFAC = 1 

NPBDC = 3 NSPFAC = 2 

ITIME Time variable. Input 0 = time independent. 

This is used for transient computation and is not 
used in steady state computations. 



Corner Node Location Cards 

Corner node number (n) 

x (n) x coordinate of node(n) 

y(n) y coordinate of node(n) 

Notes There should be NNPC number of these cards. Use 
appropriate dimensions, i.e. feet, meters which will be carried 
throughout the remaining input variables. 

NPMIS Enter a "1" or "0" where 

1 - additional nodes are to be generated after this 
node. 

0 - no nodes will be generated after this node. 
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Base Reference Elements Identification Cards 



NOD 


For each base element (N), give a counter-clockwise 
numbering for each node in the element. 




3 # I3 ' 20 

11 2 " 19 

.1 12 , 18 



NMIS 


Total number of elements above the base element before 
the next base element. 


X K < N ) 


Hydraulic conductivity of the elements defined by the 
base element in the x direction. 


YK(N) 


Hydraulic conductivity in the y direction. 


ALP (N) 


Angle of possible stratification or plane of soil 
matrix with the horizontal. 

y 

^ 1 



Note: There should be NELEMC number of these cards provided. 

Neuman Boundary Condition Cards 



NSTART 


The beginning or lowest corner node number of each 
Neuman Boundary defined. 


BVAL 


The Neuman Boundary value that defines the known 
potential height or water level measured from the 
datum reference line of the coordinate system. 


NBSAME 


Total number of boundary faces having same Neuman 
boundary condition. Each face is a line made by two 

corner nodes of the same Neuman boundary. 
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NINC 



The number of nodal number increases or numerical 
difference to the next corner node from node NSTART. 
Do not count mid-nodes only corner nodes. 



Note: There should be NF'BOC cards provided for this set. 

Dirichlet Boundary Condition Cards 



NSTAR 


The beginning or lowest corner node number which 
defines each Dirichlet boundary. 


DBVAL 


The Dirichlet boundary value that defines the potential 
or water level height measured from the data reference 
elevation of the coordinate system. 


NDSAME 


Total number of boundary faces of this Dirichlet 
boundary condition. Each face is counted between 

corner nodes defining the boundary or number of 
el ements. 


NDINC 


The number of nodal number increases or numerical 
difference to the next corner node from node NSTAR. 


NMID 


The number of nodal number increases from the NSTAR to 
the next mid-node. 



Note: There should be NPBDC cards provided for this set. 

Free Surface Nodal Cards (Maximum 16 pairs per card) 



NPFSA 


Free surface corner node number. 


NPBOT 


Node number of the bottom of the movable section of 
the grid that defines the limits of moving nodes for 
NPFSA. 



Note: There should be NF'FS pairs of node numbers given for a 

maximum of 16 pairs per card. Use as many cards as required. The 
bottom node (NF'BOT) defines the bottom of a seepage face 
correspondi ng to the grid geometry and corner node defining the 
free surface. The program uses this band of nodes and elements to 
converge on the free surface. 
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Mid-point Free Surface Nodal Cards. 



NPFSM 


The mid-point nodal numbers between corner nodes 
defining the free surface (NPFS - 1). 



Seepage Face Cards 



NSPC 


Beginning node number of the defined seepage face. 


NSSAME 


Number of seepage faces defined by elements for this 
seepage face NSPC. 


NS INC 


Node numbering increase between nodes defined by 
NSSAME. 



Note: There should be NSPFAC cards provided -for this set. 

Statement and End Card 

End of Problem Statement should be less than 72 characters and 



end with a 


peri od . 


NDIMEN 


Input a "3" to end the problem. 


ERR 


Specify any error as for FERR. 
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UNCONFINED FLOW EXAMPLE PROBLEM 




Title and Type Problem Card 



TITLE - SILTY SAND DAM 

DIMEN - 2 (two dimensional) 

FERR - 0.01 (error criterion is 0.01) 

Control Cards 



NMPC - 

NELEMC- 
NPBOC - 
NPFS - 

ITGIV - 
NPBDC - 



NSPFAC- 
ITIME - 



12 (corner nodes are 18, 52 35, 45, 52, 62, 69, 79, 86 and 
96) 

5 ( -five base elements - shaded) 

0 ( no Neuman Boundaries) 

6 ( 6 nodal points on the -free surface - 11, 28, 45, 62, 

79 and 96) 

50 ( maximum number o-f iterations) 

2 (There are two -faces with Dirichlet boundary conditions, 
the le-ft side and the right side of the mesh which 
forms the dam) 

1 (There is one seepage face - on the right side where 

flow is to be considered.) 

0 (Time independent for steady state conditions) 
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Corner Node Location Cards 



1 , 0,0, 1 

Where, the first node is number 1 having an >: coordinate of 0 and 
a y coordinate of 0. Generation of additional nodes will follow. 

11 , 0 , 110,0 

Where, node 11 is at the end of the first column at 
coordinates 0 and 110 and no generation is required for other 
nodes in this column. 



18 , 200 , 0, 1 
28 , 200 , 9? , 0 

35 . 400 . 0 . 1 

45 . 400 . 88. 0 

52 . 600 . 0 . 1 

62 . 600 . 77.0 

69 . 800 . 0 . 1 

79 . 800 . 66 . 0 

86 . 1000 . 0 . 1 

96 . 1 000 . 55 . 0 



Base Reference Elements Identification Cards 



1 , 12, 18, 19,20, 13,3,2,4,0. 002 , 0 . 002 , 0 



Counter clockwise numbering of nodes for each element and the 
permeability, 0.002 ft/day for this base element and 4 elements 
above the base element. 



18 . 29 . 35 . 36 . 37 . 30 . 20 . 19 . 4 . 0 . 002 , 

35 . 46 . 52 . 53 . 54 . 47 . 37 . 36 . 4 . 0 . 002 , 

52 . 63 . 69 . 70 . 71 . 64 . 54 . 53 . 4 . 0 . 002 , 

69 . 80 . 86 . 87 . 88 . 81 . 71 . 70 . 4 . 0 . 002 , 



0 . 002 , 0 
0 . 002 , 0 
0 . 002 , 0 
0 . 002 , 0 



Neuman Boundary Condition Cards 

None specified, so no cards will be provided 
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Dinchlet Boundary Condition Cards 

1.110.5.2.1 

36.55.5.2. 1 

The -first node o-f the first Dirichlet Boundary is node 1 and the 
DBVAL value is 110 feet. There are 5 elements on this face that 
makes this Dirichlet boundary condition, NDSAME = 5. NDINC equals 
2 for the nodal number increase for the consecutive elements, as 
in 3 - 1, 5-3, 7-5, etc. NMID equals 1 for the node numbering 

increase between successive nodes. 

Free Surface Nodal Cards 

11, 1,28, 18,45,35,62,52,79,69,96,86 

Here the total mesh is allowed to move, where the user may fix the 
bottom nodes at any level in the mesh. 

Mid-point Free Surface Nodal Cards 

17,34,51,68,85 

Seepage Face Cards 

86 , 2,2 

The seepage is prescribed to start at node 86, NSFC = 86, and 
there are 2 elements for this seepage face. NSINC is equal to 2 
as there is a two nodal increase between element corner nodes, as 
is 88 — 86 and 9 0 — 88. 

Statement and End Card 

END OF SILTY SAND DAM. 

3 , . 00 1 
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SILTY SAND DAM . 

2 , 0.01 

12.5.0. 6.50.2.1.0 

1 . 0 . 0. 1 

11 , 0 , 110,0 

18,200,0,1 

28.200.99.0 

15 . 400.0. 1 

45 . 400.88 . 0 
>2,600,0, 1 

62.600.77.0 

69.800.0. 1 

79.800.66.0 

86 . 1000 . 0 . 1 

96.1000.55.0 

1,12, 18, 19, 20, 13, 3, 2, A, 0.002, 0.002,0 
18 , 29 , 35 , 36 , 3 7 , 30 , 20 , 19 , A , 0 . 002 , 0 . 002 , 
35 , A6 , 52 , 5 3 , 5 A , A7 , 37 , 36 , A , 0 . 002 , 0 . 002 , 
52, 63, 69, 70, 71, 6A,5A, 53, A, 0.002, 0.002, 
59, 80, 86, 87, 88, 81, 71, 70, A, 0.002, 0.002, 

1.110.5.2.1 

36.55.5.2.1 

11, 1,28, 18, A5, 35, 62, 5 2, 79, 69, 96, 86 

17, 3A, 51, 68, 85 

36,2,2 

2ND OF SILTY SAND DAM. 

3,0.01 



0 

0 

0 

0 



SILTY 5 A N 0 OAH 
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47 


64 


56 


59 


46 


39 


36 


. 20000E -02 


. 2 OOOOE * 02 


0 


13 


39 


46 


59 


57 


69 


49 


4 1 


40 


. 20000E-02 


. 2 OOOOE * 02 


0 


14 


4 1 


49 


69 


59 


60 


60 


43 


42 


. 20000E * 02 


. 20000E-02 


0 


1 5 


43 


60 


90 


61 


62 


6 1 


46 


44 


. 20000E * 02 


. 2 OOOOE * 02 


0 


1 6 


62 


63 


69 


70 


7 1 


94 


64 


63 


. 20000E * 02 


. 2 OOOO E * 02 


0 


1 7 


64 


94 


7 1 


72 


73 


66 


59 


66 


. 2 OOOOE * 0 2 


. 2 OOOO E * 02 


0 


19 


69 


95 


73 


74 


79 


96 


56 


57 


. 20000E * 02 


. 20000E *02 


0 


1 9 


66 


66 


75 


79 


77 


C 7 


90 


69 


. 20000E * 02 


. 20000E *02 


0 


20 


90 


67 


77 


76 


79 


99 


92 


61 


. 20000E * 02 


. 20000E -02 


o 


2 1 












6 1 


7 1 


70 


. 20000E -02 


. 20000E -02 


0 


32 


71 


6 1 


99 


66 


60 


62 


73 


72 


. 20000E -02 


. 20000E -02 


0 


23 


73 


62 


90 


9 1 


92 


63 


75 


74 


. 2OOO0E - 02 


. 20000E -02 


0 


24 


76 


63 


92 


93 


64 


64 


77 


7 6 


. 20000E - 02 


. 2 OOOOE - 02 


0 


26 


77 


64 


94 


95 


99 


66 


79 


76 


. 20000E -02 


. 20000E -02 


0 



OIRICHLET BOUNDARY CONOITION OATA 



OES 



RESPECTIVE OIRICHLET BOUNOARY CONDITIONS 



1 

3 



2 


110. OOOO 


110. OOOO 


4 


110 . oooo 


110. OOOO 


6 


1 1 o . oooo 


1 io . oooo 


6 


110 . oooo 


110 . oooo 


1 0 


110 . oooo 


110 . oooo 


1 O 


110. oooo 


110. oooo 


67 


55 . OOOO 


56 . OOOO 


99 


55 . OOOO 


56 . OOOO 


9 1 


55 . OOOO 


65 . OOOO 


93 


66 . OOOO 


66 . OOOO 


66 


56 . OOOO 


56 . OOOO 


66 


56 . OOOO 


66 . OOOO 



) A L PO I NTS ON FREE SURFACE 

H 17 26 34 46 51 62 66 76 66 66 O 

SEEPAGE FACE NOOAL POINT OATA 
INSEOUTIYE NODE NUMBERS ON THE SEEPAGE FACE 



I' 66 

I- 60 



A- 18 



JOE * 


62 


X 


■ 


600 . 0000 


Y « 


6 1 .77 66 


JDEi 


45 


X 


■ 


400 . OOOO 


Y « 


92 . 3734 


JOE ■ 


26 


X 


■ 


200 . 0000 


Y « 


101 .7151 


JOE * 


1 1 


X ■ 0.0000 

CONVERGENCE CHECK 


Y * 

ERROR ■ 


110. OOOO 

. 1 7 1 OE - 02 



IDE 


X COORDINATE 


Y COORO 1 NATE 


V- POTENT I AL 


X - VELOCITY 


Y- VELOC I TY 


TOTAL VELD 


1 


0 . OOOO 


O . OOOO 


110. OOOO 


• . 63 1 6E 


- 0 4 


-.11 2 5 E -06 


. 63 1 6E 


- O 4 


2 


0 . OOOO 


1 1 . OOOO 


110. OOOO 


- . 83 1 5E 


- 04 


- . 1 567E - 1 5 


. 83 1 SE 


-04 


3 


O . OOOO 


22 . OOOO 


110. OOOO 


- . 6 3 0 6 E 


- 04 


. 625 1 E - 07 


. 8 3 0 6 E 


- 04 


4 


0 . OOOO 


33 . OOOO 


110. OOOO 


- . 6 2 9 6 E 


- 04 


0 . 


. 629 6E 


- O 4 


S 


0 . OOOO 


44 . OOOO 


110. OOOO 


- . 8279E 


- 04 


. 1 7S9E - 06 


. 6279E 


- 04 


« 


0 . OOOO 


55 . OOOO 


110. OOOO 


- . 6 2 5 5 E 


- 04 


0 . 


. 6255E 


- 04 


7 


0 . OOOO 


66 . OOOO 


110. OOOO 


- . 6 2 27 E 


- 04 


. 29 20E - 06 


. 6 2 2 7 E 


- 04 


« 


0 . OOOO 


77 . OOOO 


110. OOOO 


- . 6 1 9 1 E 


-04 


0 . 


. 6 1 9 1 E 


- 04 


« 


0 . OOOO 


66 . OOOO 


110. OOOO 


• . 6 1 49E 


- 04 


. 4 49 9E -06 


. 6 1 4 S E 


- 04 


10 


0 . OOOO 


9 9 . OOOO 


110. OOOO 


- . 609 9 E 


- 04 


. 1056E - 1 5 


. 6 0 9 9 E 


- 04 


1 1 


0 . OOOO 


110. OOOO 


110. OOOO 


- . 6 04 5 E 


- 04 


. 7677E -09 


. 804SE 


- 04 


12 


1 00 . OOOO 


0 . OOOO 


105 . 7709 


- . 6 600 E 


-04 


- . 2 139E -07 


. 6 9 00 E 


- O 4 


13 


1 OO . OOOO 


21.1715 


105.7743 


- . 6 59 5E 


- 04 


. 5 5 7 4 E -06 


. 6 5 9 5 E 


- 04 


1 4 


1 00 . OOOO 


42.3430 


105.7846 


- . 6562E 


-04 


. 1 1 3 1 E - 05 


. 6563E 


- O 4 


1 5 


1 OO . OOOO 


63.5145 


105.8024 


- . 6 5 9 O E 


- 04 


. 1 7 3 9 E -05 


. 6 5 6 2 E 


- 04 


1 S 


1 00 . OOOO 


64.6660 


105 . 6283 


- . 8 530E 


- 04 


. 2 4 0 3 E - 05 


. 6 5 3 3 E 


- 04 


1 7 


1 OO . OOOO 


105.6575 


105.6636 


• . 8 49 1 E 


-04 


. 3 2 2 3 E -05 


. 6 4 9 7 E 


- 04 


1 6 


200 . OOOO 


O . OOOO 


101. 4005 


• . 8 895E 


-04 


. 1 1 66E - 07 


. 669SE 


-04 


1 » 


200 . OOOO 


10. 1716 


1 O 1 .4014 


- . 8 895 E 


-04 


. 3 6 7 3 E - 06 


. 6695E 


-04 


20 


200 . OOOO 


20 . 3430 


1 O 1 . 4044 


- . 8 695E 


-04 


. 8 1 67E - 06 


. 6 695E 


-04 


2 1 


200 . OOOO 


30.5145 


101.4092 


- . 8 895E 


- 04 


. 1 1 5 7 E - 05 


. 66 9 6E 


-04 


22 


200 . OOOO 


40.6660 


101.4162 


- . 8 895E 


- O 4 


. 1 6 3 0 E - 05 


. 6696E 


-04 


23 


200 . OOOO 


50 . 6575 


101. 4246 


• . 8 69 6 E 


- 04 


. 1 9 0 9 E - 05 


. 6696E 


-04 


24 


200 . OOOO 


61 .0291 


101 .4356 


- . 6 6 9 6 E 


-04 


. 243 1 E -05 


6 6 9 9 E 


-04 


2 S 


200 . OOOO 


71 .2006 


101 .4460 


- . 669 6E 


- O 4 


. 2 6 2 2 E -05 


. 6 9 0 2 E 


- 04 


26 


200 . OOOO 


61 .3721 


101 .4622 


• . 6 6 9 9 E 


- 0 4 


. 3 2 0 6 E -05 


. 6 9 0 5 E 


-04 


27 


200 . OOOO 


9 1 .543 6 


1 O 1 .4760 


- . 6904E 


- 04 


. 3255E - 05 


. 69 1 OE 


-04 


29 


200 . OOOO 


1 O 1 .7151 


101 .4954 


- . 6 904 E 


-04 


. 39 64E -05 


. 69 1 3E 


-04 


29 


3 00 . OOOO 


O . OOOO 


96.6349 


• . 9356E 


- 04 


. 67 9 9 E - 09 


. 93 56E 


-04 


30 


300 . OOOO 


1 9 . 4066 


96.6392 


- . 9354E 


- 04 


. 673SE - 06 


. 93 5SE 


- 04 


3 1 


3 00 . OOOO 


36.6177 


96.6516 


• . 9 3 5 0 E 


-04 


. 1 74SE - 05 


. 93S2E 


-04 


32 


300 . OOOO 


66 . 226S 


96.6730 


* . 9 3 4 3 E 


-04 


.26 1 1 E - 05 


. 93 4 6E 


-04 


33 


300 . OOOO 


77 . 6354 


96.9026 


- . 9 3 3 2 E 


- 04 


. 3 4 6 6 E • 05 


. 9 3 39 E 


- 04 


34 


300 . OOOO 


97.0442 


9 6 .9406 


* . 9 3 1 6E 


- 04 


. 4 33 1 E - 05 


. 9 3 2 6 E 


- 04 


3 5 


400 . OOOO 


O . OOOO 


92.0447 


- . 9 60 1 E 


- 04 


- . 4037E - 06 


. 9 60 1 E 


-04 


3 K 


4 00 . OOOO 


9.2373 


92.0459 


- . 9 80 1 E 


- 04 


. 492 1 E -06 


. 9 60 1 E 


-04 


37 


400 . OOOO 


184747 


92. 0492 


- . 9 6 00 E 


- 04 


. 9766E • 06 


. 9 6 OO E 


-04 


3 ( 


400 . OOOO 


27.7120 


92 . 0 S 4 9 


* . 9796E 


- 04 


. 1 4 7SE - 06 


. 979 J E 


-04 


39 


400 . OOOO 


36.9494 


92.0629 


* . 9 79 SE 


- 04 


. 1 9S7E - 05 


. 9797E 


-04 


40 


400 . OOOO 


4 6. 1 667 


92.0731 


- . 9 792 E 


- 04 


. 2 4 5 4 E - 05 


. 9796E 


-04 


4 1 


400 . OOOO 


65 . 4240 


92.0655 


- . 9766E 


- 04 


. 29 33E • 05 


. 979 2E 


- 04 


4 2 


400 . OOOO 


64.6614 


92.1 002 


- . 9 763E 


- 04 


. 3 4 2 3 E -05 


. 9 7 6 9 E 


- 04 


43 


400 . OOOO 


73.6967 


•2.1172 


• . 9 7 7 6 E 


- 04 


. 3 908 E - 05 


. 9 7 A 6E 


-04 


44 


400 . OOOO 


63. 1361 


92 . 1362 


• . 9 772E 


-04 


. 4369E -05 


. 9 7 6 2 E 


-04 


45 


400 . OOOO 


92.3734 


92 . 1676 


• . 9795E 


- 04 


. 4874E-05 


. 9777E 


-04 


4 9 


500 . OOOO 


0 . OOOO 


66. 9674 


- . 1 043E 


- 03 


. 2 1 4 6 E - 06 


. 1 04 3E 


- 03 


47 


600 . OOOO 


17.4160 


66.9922 


* . 1 043E 


-03 


. 1 1 1 3E - OS 


. 1 04 3E 


-03 


49 


600 . OOOO 


34.8300 


67.0066 


* . 1 043E 


- 03 


. 2223E - OS 


. 1 04 3E 


- 03 


49 


600 . OOOO 


52.2451 


87.0306 


- . 1042E 


-03 


. 3328E -05 


. 1042E 


-03 


SO 


600 . OOOO 


6 9 . 6 60 1 


67.0642 


• . 1 040E 


- 03 


. 44 20E - 06 


. 1 04 1 E 


-03 


S 1 


600 . OOOO 


87.0751 


67 . 1073 


• . 1 039E 


- 0 3 


. 5 6 0 6 E - 08 


. 1 040E 


* 03 


S 2 


600 . OOOO 


0 . OOOO 


6 1.6116 


• . 1 1 08E 


-03 


- . 6594E-08 


. 1 1 0 6E 


-03 


S3 


800 . OOOO 


6 . 1777 


61 . 6129 


• . 1 101E 


-03 


. 6378E *06 


. 1 1 06E 


-03 



PRES . 



S 4 


600 . OOOO 


1 6 . 3 6 S 4 


61.6166 


- . 1 1 06E - 03 


. 1 270E - 06 


. 1 1 O 6 E -03 


15 


600 . OOOO 


24 . 63 30 


61.6234 


- . 1 1 08E -03 


. 1 9 1 OE - 06 


. 1 1 08E - 03 


St 


600 . OOOO 


32.7107 


61.6325 


- . 1 1 05E -03 


. 2S39E -06 


. 1 106E-03 


S7 


600 . OOOO 


40.6664 


61.6441 


- . 1 1 06E -03 


. 3 1 7 1 E - 05 


. 1 1 O 5 E -03 


St 


600 . OOOO 


49.0691 


81.6664 


- . 1 1 04E -03 


. 3 6 0 2 E - 06 


. 1 1 OSE - 03 


1» 


600 . OOOO 


57.2436 


61 .8751 


- . 1 1 03E - 03 


. 4406E - 06 


. 1 1 04E - 03 


10 


600 . OOOO 


66.4214 


61.6944 


- . 1 1 02 E -03 


. SOS3E - 05 


. 1 1 03E - 03 


1 1 


800 . OOOO 


73.5991 


81.7160 


• . 11 0 1 E - 03 


. 5 6 6 9 E - 05 


. 1 1 O 3 E -03 


12 


600 . OOOO 


6 1 .77 66 


61 .7401 


-.11 OOE - 03 


. 6295E - 05 


. 1 1 0 2 E -03 


13 


700 . OOOO 


0 . OOOO 


75.8434 


- . 1 2 0 2 E -03 


. 3 2 1 7E - 09 


. 1 2 0 2 E -03 


14 


700 . OOOO 


15. 1314 


75 . 6469 


- . 1 2 0 2 E -03 


. 1 4 6 1 E - 05 


. 1 202E - 03 


IS 


700 . OOOO 


30.2626 


75.8655 


- . 1 20 1 E - 03 


. 29 1 6E -05 


. 1 20 1 E - 03 


19 


• 700.0000 


45.3942 


75.8932 


- . 1 2 OO E - 03 


. 4 3 5 6 E -05 


. 1 20 1 E -03 


17 


700 . OOOO 


60.5256 


75.9322 


- . 1 1 99E - 03 


. 5782E - 05 


. 1 200E -03 


19 


700 . OOOO 


75 . 6570 


75.9827 


- . 1 1 9 7E - 03 


. 7233E -05 


. 1 1 99E -03 


19 


600 . OOOO 


0 OOOO 


69.6922 


- . 1 30 1 E - 03 


. 1 B 1 SE - 07 


. 1 30 1 E -03 


'0 


600 . OOOO 


6 . 9537 


69.5935 


- . 1 30 1 E -03 


. 8277E -06 


. 1 30 1 E -03 


' 1 


600 . OOOO 


1 3 . 9074 


69.5960 


- . 1 3 00 E - 03 


. 1 849 E - 05 


. 1 300E - 03 


•2 


800 . OOOO 


20.6611 


69 . 6049 


- . 1 2 9 9 E -03 


. 24 6 1 E - 05 


. 1 299E -03 


2 


800 . OOOO 


27.8149 


69.6151 


- . 1 2 9 6 E -03 


. 367 1 E -05 


. 1 2 9 6E - 03 


1 


600 . OOOO 


34 . 7686 


69 . 6275 


- . 1 296E -03 


. 4 0 2 4 E -05 


. 1 297E -03 


5 


800 . OOOO 


41.7223 


99.6431 


- . 1 2 9 4 E - 03 


. 543 1 E -05 


. 1 295E -03 


9 


600 . OOOO 


4 8.6760* 


69.6606 


• . 1 292E -03 


. 54S6E -05 


. 1 293E -03 


7 


600 . OOOO 


S8 . 62 97 


69.6610 


- . 1 2 6 6 E -03 


. 7 0 6 3 E • 05 


. 1 290E - 03 


• 


800 . OOOO 


62.5634 


69 . 7030 


- . 1 2 6 6 E -03 


. 9 670E * 05 


. 1 2 6 7 E -03 


9 


800 . OOOO 


99 . 537 1 


69.7274 


- . 1 2 6 2 E -03 


. 6 6 7 0 E -05 


. 1 2 6 5 E -03 


0 


900 . OOOO 


O . OOOO 


62 . 9 656 


- . 1 4 5 9 E -03 


- . 1 1 08E -06 


. 1 4 5 9 E -03 


1 


900 . OOOO 


12.4537 


62.6925 


- . 1 4 60E - 03 


. 1 723E - 05 


. 1 4 6 O E - 03 


2 


900 . OOOO 


24.9074 


62.7129 


- . 1 460E -03 


. 3 4 7 6 E -05 


. 1 4 6 1 E - 03 


3 


900 . OOOO 


37.3811 


62.7475 


• . 1462E -03 


. 5290E - 05 


. 1 4 6 3 E -03 


4 


900 . OOOO 


49.6149 


62.7972 


- . 1 4 64 E - 03 


. 7 1 9 8E - 05 


. 1 4 9 6E -03 


5 


900 . OOOO 


62 . 2666 


62.6634 


- . 1 4 9 6E - 03 


. 9437E - 05 


. 1 4 69 E - 03 


« 


1 OOO . OOOO 


0 . OOOO 


55 . OOOO 


- . 1 6 1 5E -03 


- . 1 767E - 06 


. 1 6 1 5E - 03 


7 


1 OOO . OOOO 


5 . 5000 


55 . OOOO 


- . 1 6 1 6 E - 03 


. 1 058E - 1 5 


. 1 8 1 6E - 03 


9 


1 OOO . OOOO 


1 1 . OOOO 


55 . OOOO 


- . 1 6 1 7E - 03 


. 502 1 E - 06 


. 1 6 1 7E * 03 


9 


1 OOO . OOOO 


1 6 . 5000 


55 . OOOO 


- . 1 620E - 03 


- . 1 OS6E - 1 5 


. 1 6 20E - 03 


0 


1 OOO . OOOO 


2 2 . OOOO 


55 . OOOO 


- . 1 623E -03 


. 1 03 4 E * 05 


. 1 623E - 03 


1 


1 OOO . OOOO 


2 7 . SOOO 


55 . OOOO 


- . 1 629E -03 


0 . 


. 1 629E - 03 


2 


1 OOO . OOOO 


33 . OOOO 


55 . OOOO 


- . 1 634E -03 


. 1 627E - 05 


. 1 634 E - 03 


3 


1 OOO . OOOO 


36 . SOOO 


55.0000 


- . 1 642E -03 


O . 


. 1 B42E - 03 


4 


1 OOO . OOOO 


4 4 . OOOO 


55 . OOOO 


- . 1 649E -03 


. 231 9E -OS 


. 1 BSOE - 03 


5 


1 OOO . OOOO 


4 9 . SOOO 


55 . OOOO 


- . 1 66 1 E -03 


0 


. 1 6 6 1 E - 03 


9 


1 OOO . OOOO 


S 5 . OOOO 


S 5 . OOOO 


- . 1 696E-03 


. 3326E -OS 


. 1 6 66E - 03 



SEEPAGE FACE COMPUTATION RESULTS 



6 6 


NE s 


6 6 


AVE . 


VE L * 


1 6 1 60E - 03 


OS 3 


. 1 1 O OOE ♦ O 2 


DO > 


. 1 7779E -02 


TOTAL 


D 1 SCHARGE 


i 


. 1 777 9 E 


-02 










8 6 


NE = 


90 


A VE . 


VE L * 


1 6202E - 03 


DS > 


. 1 1 0 OOE + 02 


DO » 


. 1 7822E -02 



TOTAL DISCHARGE ■ .3SS98E-02 



END OF SILTV SAND DAM. 

DIMEN * 3 FERR • . O 1 OOO 

A- 19 



COEF . 



APPENDIX B— 1 



SAMPLE INPUT 



AND OUTPUT 



RUN ADS 2 SP = 25 FEET DEPTH = 20 FEET K = .2. DEPTH =58'. 

2,0.0001 

73. 30. 0. 30. 20. 6. 5.0 

I, -340,0,1 
5,-340,34,1 

II, -340,78.0,0 
18,-200,0,1 
22,-200,34,1 
28,-200,78.0,0 
15,-80,0, 1 
19,-80,34,1 
'45,-80,77.0,0 

52.0. 0.1 

56.0. 34.1 

52.0. 74.0.0 

59.50.0. 1 

73.50.34.1 

79.50.71.0. 0 

16.58.75.0. 1 

70.58. 75.34.1 

76.59.0. 58.0.0 

103.62.0. 0.1 

107.62.0. 34.1 

113.61.0. 58.0.0 

120.66.0. 0.1 

124.66.0. 34.1 

130.65.0. 60.0.0 

137.70.0. 0.1 

141.70.0. 34.1 

147.69.0. 60.0.0 

154.74.0. 0.1 

158.74.0. 34.1 

164.75.0. 60.0.0 

171. 78. 0. 0.1 

175.78.0. 34.1 

181.79.0. 60.0.0 

188.82.0. 0.1 

192.82.0. 34.1 

198.83.0. 58.0.0 
£05,86.0,0,1 
>09,86.0,34,1 
>15,85.0,58.0,0 
>22,90.0,0,1 
>26,90.0,34,1 
>32,89.0,60.0,0 
>39,94.0,0,1 
>43,94.0,34,1 
>49,93.0,60.0,0 

156.98.0. 0.1 

160.98.0. 34.1 

166.99.0. 60.0.0 

173.102.0. 0.1 

177.102.0. 34.1 

183.103.0. 60.0.0 

190.106.0. 0.1 

194.106.0. 34.1 

100.107.0. 58.0.0 

107.110.0. 0.1 

111.110.0. 34.1 

117.109.0. 58.0.0 



B-l 



124.114.0. 0.1 

128.114.0. 34.1 

134.113.0. 60.0.0 

141.118.0. 0.1 

145.118.0. 34.1 

151.117.0. 60.0.0 

158.122.0. 0.1 

162.122.0. 34.1 

168.123.0. 60.0.0 

175.126.0. 0.1 

179.126.0. 34.1 

185.127.0. 60.0.0 

192.130.0. 0.1 

196.130.0. 34.1 
>02, 131.0,58.0,0 
>09,134.0,0,1 
>13,134.0,34,1 
>19,133.0,58.0,0 
>26,138.0,0,1 
>30,138.0,34,1 

36. 137.0. 60.0.0 

43.142.0. 0.1 
•47,142.0,34,1 

53.141.0. 60.0.0 

60.146.0. 0.1 

64.146.0. 34.1 

70.147.0. 60.0.0 

77.150.0. 0.1 
•81,150.0,34,1 
•87,151.0,60.0,0 

94.154.0. 0.1 

98.154.0. 34.1 

04.155.0. 58.0.0 

11.156.0. 0.1 

15.156.0. 34.1 

21.156.0. 58.0.0 

,12,18,19,20,13,3,2,4,0.2,0.2,0.0 

8.29.35.36.37.30.20.19.4.0. 2.0.2.0.0 

5.46.52.53.54.47.37.36.4.0. 2.0.2.0.0 

2.63.69.70.71.64.54.53.4.0. 2.0.2.0.0 

9.80.86.87.88.81.71.70.4.0. 2.0.2.0.0 

6.97.103.104.105.98.88.87.4.0. 2.0.2.0.0 

03.114.120.121.122.115.105.104.4.0. 2.0.2.0. 

20.131.137.138.139.132.122.121.4.0. 2.0.2.0. 

37.148.154.155.156.149.139.138.4.0. 2.0.2.0. 

54.165.171.172.173.166.156.155.4.0. 2.0.2.0. 

71.182.188.189.190.183.173.172.4.0. 2.0.2.0. 

88.199.205.206.207.200.190.189.4.0. 2.0.2.0. 

05.216.222.223.224.217.207.206.4.0. 2.0.2.0. 

22.233.239.240.241.234.224.223.4.0. 2.0.2.0. 

39.250.256.257.258.251.241.240.4.0. 2.0.2.0. 

56.267.273.274.275.268.258.257.4.0. 2.0.2.0. 

73.284.290.291.292.285.275.274.4.0. 2.0.2.0. 

90.301.307.308.309.302.292.291.4.0. 2.0.2.0. 

07.318.324.325.326.319.309.308.4.0. 2.0.2.0. 

24.335.341.342.343.336.326.325.4.0. 2.0.2.0. 

41.352.358.359.360.353.343.342.4.0. 2.0.2.0. 

58.369.375.376.377.370.360.359.4.0. 2.0.2.0. 

75.386.392.393.394.387.377.376.4.0. 2.0.2.0. 

92.403.409.410.411.404.394.393.4.0. 2.0.2.0. 

B- 



,0 

,0 

,0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

-2 



409, 420, 426,427, 428, 421, All, 410, 4, 0.2, 0.2, 0.0 

426.437.443.444.445.438.428.427.4.0. 2.0.2.0.0 

443.454.460.461.462.455.445.444.4.0. 2.0.2.0.0 

460.471.477.478.479.472.462.461.4.0. 2.0.2.0.0 

477.488.494.495.496.489.479.478.4.0. 2.0.2.0.0 

494.505.511.512.513.506.496.495.4.0. 2.0.2.0.0 

1.78.0. 5.2.1 

96.58.0. 1.17.6 

198.58.0. 1.17.6 

300.58.0. 1.17.6 

402.58.0. 1.17.6 

504.58.0. 1.17.6 

11,5,28,22,45,39,62,56,79,73,96,90,113,107,130,124 
147,141,164,158,181,175,198,192,215,209, 232 , 226 , 249 , 243 , 266 , 260 
283,277,300,294,317,311,334,328,351,345,368,362,385,329,402,396 
419,413,436,430,453,447,470,464,487,481,504,498 
17,34,51,68,85,102,119,136,153,170,187,204,221,238,255,272 
289 , 306 ,323,340,357,374,391, 408 , 425 , 442 , 459 , 476 , 493 
96,1,17,198,1,17,300,1,17,402,1,17,504,1,17 
SP-25..K = .2..D = 20 FT... END OF PROBLEM. 

3.0. 0001 



_ 



B-3 



1 



RUN ADS 2 



5 p • 25 Feet 



DEPTH = 20 FEET K * 2 



DEPTH * 56' 



NDIMEN • 



NNPC • 63 



N E L EMC ■ 30 



NPF5 * 30 



ITGIY • 20 



NP5DC 



N5PFAC ■ 



I T I HE 



GENERATED NDOAL POINT DATA AND (X)-(Y) CDDRDINATES 



NODE 


1 


X * 


- 340 . OOOO 


Y B 


0 . OOOO 


NDDE 


3 


X ■ 


- 340 . OOOO 


Y b 


1 7 . OOOO 


NODE 


5 


X ■ 


- 340.0000 


Yb 


34 . OOOO 


NDDE 


7 


X ■ 


- 340 . OOOO 


Y B 


48 . 6557 


NDDE 


6 


X ■ 


•340.0000 


Yb 


63 . 3333 


NDDE 


1 1 


X ■ 


- 340.0000 


Yb 


78 . OOOO 


NDDE 


1 6 


X ■ 


• 200 . OOOO 


Y« 


0 . OOOO 


NDDE 


20 


X > 


- 200 . OOOO 


Yb 


1 7 .0000 


NDDE 


22 


X ■ 


- 200 . OOOO 


Yb 


34 .OOOO 


NDDE 


24 


X ■ 


- 200 . OOOO 


Yb 


46.6667 


NDDE 


26 


x» 


-200.0000 


Yb 


63 . 3333 


NDDE 


25 


X ■ 


-200.0000 


Y b 


76 . OOOO 


NDDE 


35 


X ■ 


•60.0000 


Y b 


0 . OOOO 


NDDE 


37 


X « 


-60 OOOO 


Y b 


1 7 . OOOO 


NDDE 


36 


x« 


-60 OOOO 


Y b 


34 . OOOO 


NDDE 


4 1 


X ■ 


• 60 . OOOO 


Y b 


46 . 3333 


NDDE 


43 


X » 


-60.0000 


Y b 


62 . 6667 


NDDE 


45 


X • 


-60.0000 


Yb 


77 . OOOO 


NDDE 


52 


X • 


O . OOOO 


Yb 


0 . OOOO 


NDDE 


54 


X* 


0 . OOOO 


Yb 


1 7 . OOOO 


NODE 


56 


X* 


O . OOOO 


Yb 


34 . OOOO 


NDDE 


56 


X* 


O . OOOO 


Y b 


47 . 3333 


NDDE 


60 


X* 


O . OOOO 


Y « 


80.6667 


NDDE 


62 


X • 


O . OOOO 


Y b 


74 . OOOO 


NDDE 


89 


X 1 


50 . OOOO 


Y b 


0 . OOOO 


NDDE 


7 1 


X • 


50 . OOOO 


Y b 


1 7 . OOOO 


NDDE 


73 


X S 


50 . OOOO 


Y B 


34 .0000 


NDDE 


75 


X B 


SO . OOOO 


Y B 


46 . 3333 


NDDE 


77 


X B 


SO . OOOO 


Y B 


56.6667 


NDDE 


79 


X B 


50 . OOOO 


Y B 


71 .OOOO 


NDDE 


65 


X B 


56 . 7500 


Y B 


O . OOOO 


NDDE 


66 


X B 


56 . 7500 


Y B 


1 7 . OOOO 


NDDE 


90 


X B 


66 . 7500 


Y B 


34 . OOOO 


NDDE 


92 


X B 


55 . 5333 


Y B 


42 . OOOO 


NDDE 


94 


X > 


55.9167 


Y B 


SO . OOOO 


NDDE 


6 6 


X B 


56 . OOOO 


Y B 


56 . OOOO 


NDDE 


1 03 


x> 


62 . OOOO 


Y B 


O . OOOO 


NDDE 


1 OS 


X B 


62 . OOOO 


Y B 


1 7 . OOOO 


NDDE 


1 07 


X B 


52 . OOOO 


Y B 


34 . OOOO 


NDDE 


1 06 


X B 


6 1 . 5667 


Y B 


42 . OOOO 


NDDE 


1 1 1 


X B 


61 . 3333 


Y B 


SO . OOOO 


NDDE 


1 13 


X B 


6 1 . OOOO 


Y B 


56 . OOOO 


NDDE 


1 20 


X B 


66 . OOOO 


Y B 


O . OOOO 


NDDE 


122 


X B 


66 . OOOO 


Y B 


1 7 . OOOO 


NDDE 


124 


X B 


66 . OOOO 


Y B 


34 . OOOO 


NDDE 


126 


X B 


66 . 6667 


Y B 


42 . 6667 


NDDE 


126 


X B 


65 .3333 


Y B 


5 1 . 3333 


NDDE 


1 30 


X B 


65 .OOOO 


Y B 


60 . OOOO 


NDDE 


1 37 


X B 


70 . OOOO 


Y B 


O . OOOO 


NDDE 


136 


X B 


70 . OOOO 


Y B 


1 7 . OOOO 


NDDE 


1 4 1 


X B 


70 . OOOO 


Y B 


34 . OOOO 


NDDE 


143 


X B 


66.6667 


Y B 


42.6667 


NDDE 


146 


x« 


69 . 3333 


Y B 


5 1 . 3333 


NDDE 


1 47 


XB 


66 . OOOO 


Y B 


60 . OOOO 


NDDE 


154 


x« 


74 . OOOO 


Y B 


0 . OOOO 


NDDE 


155 


XB 


74 . OOOO 


Yb 


1 7 . OOOO 


NDDE 


166 


XB 


74 . OOOO 


yb 


34 . OOOO 


NDDE 


1 60 


X * 


74 . 3333 


Y B 


42.6667 


NDDE 


1 62 


X B 


74 .6667 


Y B 


5 1 . 3333 


NDDE 


1 64 


X B 


76 .OOOO 


Yb 


60 . OOOO 


NDDE 


171 


XB 


76 .OOOO 


Y B 


0 . OOOO 


NDDE 


1 73 


XB 


76 . OOOO 


Yb 


17 .OOOO 


NDDE 


176 


XB 


76 . OOOO 


Y B 


34 . OOOO 


NODE 


1 77 


X* 


76 .3333 


Y B 


42 . 5657 


NODE 


176 


X* 


76 . 6667 


Y B 


5 1 . 3333 


NODE 


1 6 1 


X B 


76 .OOOO 


Yb 


60 . OOOO 


NDDE 


1 46 


XB 


62 . OOOO 


Yb 


0 . OOOO 


NDDE 


160 


XB 


62 .OOOO 


Yb 


1 7 . OOOO 


NODE 


192 


XB 


62 . OOOO 


Y b 


34 . OOOO 


NODE 


194 


X B 


62 . 3333 


Yb 


42.0000 


NODE 


16 6 


X B 


62 . 6667 


Yb 


50 . OOOO 


NODE 


166 


XB 


63 . OOOO 


Yb 


56 . OOOO 


NODE 


205 


X B 


66 .OOOO 


Y« 


0 . OOOO 


NODE 


207 


X B 


66 . OOOO 


Y b 


1 7 . OOOO 


NDDE 


209 


X B 


66 .OOOO 


Y B 


34 . OOOO 


NDDE 


21 1 


X B 


55 . 6667 


Yb 


42 . OOOO 


NODE 


213 


X B 


65 .3333 


Y B 


50 . OOOO 


NDDE 


2 15 


XB 


55 . OOOO 


Yb 


55 . OOOO 


NDDE 


222 


X B 


60 . OOOO 


Yb 


0 . OOOO 


NDDE 


224 


X a 


60 . OOOO 


Yb 


1 7 . OOOO 


NDDE 


226 


XB 


90 . OOOO 


Yb 


34 . OOOO 


NDDE 


226 


XB 


69 . 5657 


Yb 


42 . 6 667 


NDDE 


2 30 


XB 


56 . 3333 


Yb 


5 1 . 3333 


NDDE 


232 


XB 


66 . OOOO 


Y B 


60 . OOOO 


NODE 


236 


XB 


64 . OOOO 


Y b 


O . OOOO 


NDDE 


24 1 


X B 


94 . OOOO 


Y B 


1 7 . OOOO 


NDDE 


243 


XB 


64 .OOOO 


Yb 


34 . OOOO 


NDDE 


246 


XB 


63 . 6667 


Yb 


42 . 6657 


NDDE 


247 


XB 


63 . 3333 


Y b 


5 1 . 3333 


NDDE 


246 


X B 


63 . OOOO 


Y B 


60 . OOOO 


NDDE 


256 


X B 


65 .OOOO 


Y * 


0 . OOOO 


NDDE 


256 


X B 


96 . OOOO 


Y B 


1 7 . OOOO 


NDDE 


280 


X B 


66 . OOOO 


Y B 


34 . OOOO 


NDDE 


252 


X B 


96 . 3333 


Y B 


42 . 6657 


NDDE 


264 


XB 


66.6657 


Y B 


5 1 . 3333 


NDDE 


265 


X B 


66 . OOOO 


Y B 


60 . OOOO 


NDDE 


273 


X > 


102 . OOOO 


Y B 


0 . OOOO 


NDDE 


275 


X B 


102 . OOOO 


Y B 


1 7 . OOOO 


NDDE 


277 


X B 


102 . OOOO 


Y B 


34 . OOOO 


NDDE 


278 


X B 


102 . 3333 


Y B 


42 . 5657 


NDDE 


26 1 


XB 


102 . 5657 


Y B 


5 1 . 3333 


NDDE 


263 


X > 


103 .OOOO 


Y B 


60 . OOOO 


NDDE 


260 


X B 


106 . OOOO 


Y B 


0 . OOOO 


NDDE 


292 


X B 


106 . OOOO 


Y ■ 


1 7 . OOOO 


NDDE 


294 


X B 


106 . OOOO 


Y B 


34 . OOOO 


NDDE 


265 


X B 


106 . 3333 


Y B 


42 . OOOO 


NDDE 


296 


X • 


105 . 6657 


Y B 


50 . OOOO 


NDDE 


300 


X B 


1 07 . OOOO 


Y B 


58 . OOOO 


NDDE 


307 


X B 


110. OOOO 


Y B 


O . OOOO 


NDDE 


306 


X B 


110. OOOO 


Y B 


1 7 . OOOO 


NDDE 


3 1 1 


X B 


110. OOOO 


Y B 


34 . OOOO 


NDDE 


3 1 3 


X B 


1 06 . 6667 


Yb 


42 . OOOO 


NDDE 


315 


X B 


1 06 . 3333 


Y B 


SO . OOOO 


NDDE 


3 1 7 


X B 


1 09 . OOOO 


Y B 


56 . OOOO 


NDDE 


324 


X B 


114. OOOO 


Y B 


0 . OOOO 


NDDE 


326 


X B 


1 1 4 .OOOO 


Y B 


1 7 . OOOO 


NDDE 


325 


X> 


114. OOOO 


Y B 


34 . OOOO 


NDDE 


330 


X* 


1 13. 6667 


Y B 


42 . 6667 


NDDE 


332 


X B 


1 13. 3333 


Y B 


5 1 . 3333 















B-4 



NODE 


345 


X > 


1 1 6 OOOO 


V s 


34 OOOO 


NODE 


347 


X > 


1176667 


Y a 


42 6667 


NODE 


341 


X > 


1 1 7 .3333 


Y a 


5 1 .3333 


N 00 E 


36 1 


X « 


117. 0000 


Y a 


60 .0000 


NODE 


365 


X « 


122. OOOO 


Y a 


0 . OOOO 


NODE 


360 


X « 


122. OOOO 


Y * 


17 .0000 


NOOE 


362 


X > 


1 22 . OOOO 


Y * 


34 . OOOO 


NODE 


364 


X a 


122. 3333 


Y * 


42 . 6657 


NODE 


366 


X a 


122. 6667 


Y a 


6 1 . 3333 


NODE 


366 


X a 


123. OOOO 


Y a 


60 . OOOO 


NODE 


376 


X « 


126.0000 


Y a 


0 . OOOO 


NDOE 


377 


X« 


126 .0000 


Y a 


1 7 . OOOO 


NODE 


376 


X a 


126.0000 


Y a 


34 . OOOO 


NDOE 


351 


X a 


126. 3333 


Y « 


42 . 6667 


NODE 


353 


X a 


126. 6667 


Y a 


6 1 .3333 


NODE 


355 


X a 


127. OOOO 


Y « 


50 . OOOO 


NDOE 


352 


X a 


1 30 . OOOO 


Y a 


0 . OOOO 


NODE 


354 


X« 


130.0000 


Y a 


1 7 . OOOO 


NDOE 


355 


X a 


1 30.0000 


Y a 


34 . OOOO 


NODE 


366 


X a 


130.3333 


Y a 


42 . OOOO 


NOOE 


400 


X a 


130.6667 


Y a 


50 . OOOO 


NOOE 


402 


X« 


1 3 1 .0000 


Y a 


58 . OOOO 


NOOE 


406 


X a 


134. OOOO 


Y a 


0 . OOOO 


NOOE 


41 1 


X a 


1 34 . OOOO 


Y a 


1 7 . OOOO 


NDOE 


413 


X a 


134 . OOOO 


Y a 


34 . OOOO 


NODE 


4 1 5 


X a 


133 . 6667 


Y a 


42 . OOOO 


NODE 


4 17 


X a 


133 . 3333 


Y a 


60 . OOOO 


NODE 


4 1 6 


X a 


133. OOOO 


Y a 


66.0000 


NODE 


426 


X a 


138. OOOO 


Y a 


0 . OOOO 


NOOE 


426 


X a 


138. OOOO 


Y a 


1 7 . OOOO 


NOOE 


430 


X > 


138. OOOO 


Y a 


34 . OOOO 


NOOE 


432 


X > 


137.8667 


Y * 


42 . 6667 


NODE 


434 


X a 


1 37 . 3333 


Y a 


5 1 .3333 


NODE 


436 


X a 


1 37 . OOOO 


Y a 


60 . OOOO 


NODE 


443 


X a 


142 . OOOO 


Y a 


O . OOOO 


NODE 


445 


X a 


142. OOOO 


Y 5 


1 7 . OOOO 


NODE 


447 


X a 


142. OOOO 


Y 5 


34 . OOOO 


NOOE 


446 


X a 


14 1. 6667 


Y a 


42.6667 


NOOE 


45 1 


X a 


14 1 .3333 


Y a 


5 1 .3333 


NOOE 


463 


X a 


1 4 1 .OOOO 


Y a 


60 . OOOO 


NDOE 


460 


X a 


146. OOOO 


Y a 


O . OOOO 


NOOE 


462 


X a 


146. OOOO 


Y a 


1 7 . OOOO 


NOOE 


464 


X a 


146. OOOO 


Y a 


34 . OOOO 


NDOE 


466 


X a 


146.3333 


Y a 


42.6667 


NODE 


468 


X a 


146. 6667 


Y a 


6 1 .3333 


NODE 


470 


X a 


147. OOOO 


Y a 


60 . OOOO 


NODE 


477 


X a 


1 50 . OOOO 


Y a 


O . OOOO 


NODE 


475 


X a 


150. OOOO 


Y a 


1 7 . OOOO 


NODE 


46 1 


X > 


1 60 . OOOO 


Y a 


34 . OOOO 


NODE 


463 


X a 


150.3333 


Y a 


42.6667 


NODE 


466 


X a 


150.6667 


Y a 


6 1 .3333 


NODE 


467 


X > 


151. OOOO 


Y a 


60 . OOOO 


NODE 


464 


X a 


1 54 . OOOO 


Y a 


O . OOOO 


NDDE 


466 


X « 


154. OOOO 


Y a 


1 7 . OOOO 


NDDE 


466 


X * 


164. OOOO 


Y a 


34 . OOOO 


NDDE 


600 


X a 


164.3333 


Y a 


4 2 . OOOO 


NDDE 


602 


X a 


164.6667 


Y a 


50 . OOOO 


NODE 


604 


X a 


165. OOOO 


Y a 


66 . OOOO 


NDDE 


5 1 1 


X a 


166. OOOO 


Y a 


0 . OOOO 


NDDE 


613 


X a 


166. OOOO 


Y a 


1 7 . OOOO 


NODE 


6 1 6 


X a 


166. OOOO 


Y a 


34 . OOOO 


NOOE 


6 17 


X a 


156. OOOO 


Y a 


42 . OOOO 


NDDE 


616 


X a 


155. OOOO 


Y a 


60 . OOOO 


NDDE 


62 1 


X * 


156. OOOO 


Y a 


66 .OOOO 



1 NUMBER DF NODAL POINTS... S21 



NUMBER DF ELEMENTS ISO 

ELE. NO. NODAL P0INT5 



1 


1 


- 12 


1 8 


16 


20 


1 3 


3 


2 


. 2 00 00 E + 0 0 


. 2 OOOOE +00 


O 


2 


3 


13 


20 


21 


22 


1 4 


6 


4 


. 20000E+00 


. 2 OOOOE +00 


O . 


3 


6 


14 


22 


23 


24 


1 5 


7 


6 


. 200006 + 00 


. 2 OOOOE +00 


O . 


4 


7 


16 


24 


26 


26 


1 6 


6 


6 


. 2 0 000 E + OO 


. 20000E + 00 


0 . 


5 


6 


1 6 


26 


27 


26 


17 


1 1 


IO 


. 20000E+00 


. 2 OOOOE + 00 


O . 


6 


1 6 


25 


36 


36 


37 


30 


20 


1 6 


. 20000E+00 


. 2 OOOOE ♦ OO 


0. 


7 


20 


30 


37 


36 


36 


31 


22 


2 1 


. 20000E ♦ OO 


. 20000E+00 


0. 


8 


22 


31 


36 


40 


4 1 


32 


24 


23 


. 20000E ♦ 00 


. 2 OOOOE + OO 


O . 


6 


24 


32 


4 1 


42 


43 


33 


25 


26 


. 20000 E ♦ OO 


. 20000E + 00 


0 . 


io 


25 


33 


43 


44 


46 


34 


26 


27 


. 20000E + 00 


. 2 OOOOE ♦ 00 


0 . 


1 1 


35 


46 


62 


63 


54 


47 


37 


36 


. 20000E + 00 


. 2 OOOOE ♦ 0 0 


0 . 


1 2 


37 


47 










35 


36 


. 20000E + 00 


. 20OOOE + 00 


0 . 


1 3 


36 


46 


66 


67 


65 


45 


4 1 


40 


. 2000OE+ 00 


. 20000E+00 


0 . 


1 4 


4 1 


46 


66 


69 


50 


60 


43 


42 


. 200006+00 


. 20000E + 00 


0 . 


1 5 


43 


60 


60 


6 1 


62 


61 


46 


44 


. 2OO0OE +00 


. 20000E + 00 


0 . 


1 6 


62 


63 


86 


70 


7 1 


64 


54 


53 


. 20000E+00 


. 20000E + 00 


0 . 


17 


64 


64 


7 1 


72 


73 


66 


66 


56 


. 20000E+00 


. 2 OOOOE + OO 


0 . 


1 6 


66 


66 


73 


74 


76 


51 


65 


57 


. 20000E+00 


. 2 OOOOE +00 


0 . 


16 


66 


56 


76 


76 


77 


67 


60 


55 


. 2 OOOOE ♦ OO 


. 2 OOOOE + 00 


0. 


20 


60 


67 


77 


76 


76 


66 


62 


6 1 


. 2 00 00 E ♦ 00 


. 20000E + 00 


o . 


2 1 


66 


60 


66 


67 


68 


8 1 


7 1 


70 


. 20000E + 00 


. 20000E + 00 


o . 


22 


7 1 


6 1 


66 


86 


60 


82 


73 


72 


. 20000E+00 


. 20000E+00 


0. 


23 


73 


62 


90 


9 1 


92 


63 


76 


74 


. 20000E + 00 


. 20000E+00 


o . 


24 


75 


63 


92 


53 


64 


64 


77 


76 


. 2 OOOOE +00 


. 20000E + 00 


o . 


25 


77 


64 


64 


96 


96 


66 


79 


76 


. 20000E+00 


. 20000E+00 


0 . 


26 


66 


67 


1 03 


104 


105 


66 


86 


67 


. 20000E + 00 


. 20000E+00 


0 . 


27 


66 


66 


1 06 


106 


107 


66 


90 


89 


. 2 00 00 E ♦ 00 


. 2 OOOOE + 00 


o . 


28 


60 


96 


107 


108 


109 


1 00 


92 


6 1 


. 2 OOOOE ♦ OO 


. 20000E + 00 


o . 


26 


92 


ioo 


106 


1 IO 


1 1 1 


101 


64 


63 


. 2 OOOOE + 00 


. 2 OOOOE +O0 


o . 


30 
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NOOEa 11 


X a 


- 340 . OOOO 


Y a 


7ft .OOOO 
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NEW CO ORO I NATES OF THE FREE SURFACE LINE 



NOOEa 487 


X a 


1 AO .9255 


Y * 


Aft . 1 4 1 2 


NOOE a 470 


X > 


14ft. 9 2 ft 9 


Y > 


Aft . 177ft 


NOOE a 4A3 


X * 


14 1. 0 A 9 9 


Y a 


Aft . 1 S2 1 


NOOEa 438 


X > 


137.0712 


Y a 


Aft . 14 68 


NOOEa 4 1 9 


X > 


1 33 . OOOO 


Y a 


ftft . OOOO 


NOOEa 402 


X ■ 


131. OOOO 


Y a 


ft A . OOOO 


NOOEa Sftft 


X ■ 


125 . A023 


Y a 


• ft . 1S49 


NOOE a 3ftft 


X ■ 


122.931ft 


Y a 


Aft . 2272 


NOOE > 3A 1 


X > 


119.7315 


Y a 


ftft . 30 3 A 


NOOE a 334 


X a 


1 1 3 . OS93 


Y a 


Aft . 1 ftft 1 


NOOE a 3 1 7 


X a 


1 Oft . OOOO 


Y a 


5ft . OOOO 


NOOEa 300 


X a 


1 07 . OOOO 


Y « 


5ft . OOOO 


NOOE > 2B3 


X a 


102 . 9334 


Y a 


•ft . 2S72 


NOOE a 2ftft 


X * 


• ft . 9 3 • • 


Y a 


ftft . 3ft 1 4 


NOOE a 24ft 


X a 


93 .0922 


Y a 


Aft . 3ft 3 ft 


NOOEa 232 


X a 


•9 . 0943 


Y a 


ftft .329 1 


NOOE a 2 1 ft 


X a 


•ft . OOOO 


Y a 


Aft . OOOO 


NOOEa 1 ftft 


X a 


•3 . OOOO 


Y a 


5ft . OOOO 


NOOE * 1 ft 1 


X a 


7ft . 94 10 


Y a 


5ft . 4SA2 


NOOE * 1 84 


X a 


74.9499 


Y a 


5ft . 9 1 ftft 


NOOE > 147 


X a 


5ft . Oft 03 


Y a 


5ft S934 


NOOE * 1 30 


X > 


•ft . 0538 


Y > 


5ft . ftOftft 


N00E> 1 1 3 


X a 


ft 1 . OOOO 


Y > 


•ft . OOOO 


NOOE* 99 


X a 


•ft . oooo 


Y > 


•ft . oooo 


NOOEa 79 


X a 


•o . oooo 


Y > 


ftft . ftft 1 9 


NOOEa 92 


X a 


o . oooo 


Y ■ 


A 2 . SSftft 


NOOEa 4ft 


X a 


-•o . oooo 


Y > 


99 . S2 1 • 


NOOEa 2ft 


X a 


• 200 . OOOO 


Y > 


72 . 1 27A 


NOOEa 11 


X a 


• 340 . OOOO 


Y > 


7ft . OOOO 
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NEW COORDINATES OF THE FREE SURFACE LINE 



NOOE*4ft7 


X a 


150. 92ftft 


Y 


i 


5ft 


NOOEa 470 


X a 


149 9 2 ftft 


Y 


i 


5ft 


NOOEa 4A3 


X a 


141. 0999 


Y 


i 


5ft 


NOOEa 439 


X a 


137.0712 


Y 


a 


5ft 


NOOEa 4 1ft 


X a 


133 . OOOO 


Y 


a 


5ft 


NOOE *402 


X a 


131. OOOO 


Y 


a 


5ft 


NOOE a 3ftft 


X a 


1 2 A . 901 5 


Y 


a 


Aft 


NOOE*3lft 


X ■ 


122.931ft 


Y 


a 


5ft 


NOOEa 3ft 1 


X a 


lift. 7340 


Y 


a 


Aft 


NOOEa 334 


X > 


1 1 3 .OSAO 


Y 


a 


ftft 


NOOEa 317 


X > 


109 . OOOO 


Y 


a 


• ft 


NOOEa 300 


X ■ 


107 . OOOO 


Y 


a 


ftft 


NOOEa 2A3 


X ■ 


102 . ft 333 


Y ' 


a 


ftft 


NOOE»2ftft 


X « 


95 . S3ftft 


Y 


a 


• 5 


NOOEa 249 


X « 


93 . 0 ft 2 2 


Y 


a 


5ft 


NOOEa 232 


X « 


59 . 0 ft 4 3 


Y 


a 


55 


NOOEa 2 1 5 


X > 


55 .0000 


Y 


a 


5ft 


NOOEa 1 ftft 


X > 


53 . OOOO 


Y 


a 


Aft 


NOOEa 1 ftl 


X a 


7ft . 94 10 


Y 


a 


5ft 


NOOEa 1 ft4 


X a 


74 .9499 


Y 


a 


• ft 


NOOEa 147 


X a 


99 . 0 ft 04 


Y 


a 


• ft 


NOOEa 1 30 


X a 


95 . 0539 


Y 


a 


Aft 


NOOEa 1 1 3 


X > 


5 1 . OOOO 


Y 


a 


5ft 


NOOEa ft 9 


X a 


59 . OOOO 


Y 


a 


ftft 


NOOEa 7ft 


X > 


AO .OOOO 


Y 


a 


58 


NOOE* 92 


X > 


0 . OOOO 


Y 


a 


82 


NOOE* 45 


X a 


-SO . OOOO 


Y 


a 


56 


NOOE* 25 


X « 


-200 . OOOO 


Y 


a 


72 


NOOEa 11 


X a 


- 340 . OOOO 


Y 


a 


7ft 




CONVERGENCE CHECK 


ERROR 


a 




ON 0 0 E X 


COORDINATE Y COORDINATE 




1 


- 340 . OOOO 


0 . OOOO 




2 


- 340 . OOOO 


A . SOOO 




3 


-340 . OOOO 


1 7 . OOOO 




4 


- 340 . OOOO 


25 . 5000 




5 


- 340 . OOOO 


34 . OOOO 




9 


- 340 . OOOO 


4 1 . 3333 




7 


- 340 . OOOO 


45 . ft ft S 7 




ft 


- 340 . OOOO 


55 OOOO 




9 


- 340 . OOOO 


S3 . 3333 




1 O 


- 340 . OOOO 


70 . ft ft #7 




1 1 


- 340 . OOOO 


7 A . OOOO 




1 2 


- 270 . OOOO 


O . OOOO 




1 3 


- 270 . OOOO 


1 7 . OOOO 




1 4 


- 270 . OOOO 


34 . OOOO 




1 5 


- 270 . OOOO 


47 . Aft 73 




1 5 


- 270 . OOOO 


8 1 . 374ft 




1 7 


- 270 . OOOO 


75 . Oft 1 9 




1 ft 


- 200 . OOOO 


O . OOOO 




1 ft 


- 200 . OOOO 


ft . 5000 




20 


- 200 . OOOO 


1 7 . OOOO 




2 1 


- 200 . OOOO 


25 . 5000 




22 


- 200 . OOOO 


34 . OOOO 




23 


-200 . OOOO 


40. 3 ft 40 




24 


- 200 . OOOO 


45 . 7079 




25 


- 200 . OOOO 


•3.0519 





14 10 
1772 

1 ft 1 A 
1 49ft 

0000 
0000 

1 ft34 
227 1 
239 9 
204ft 
0000 
0000 
2ft7 1 
3508 
3 ft 2 2 
3250 
OOOO 
OOOO 
4955 
9 193 
• 909 
9 O ft 2 
OOOO 

0000 
57ft4 
5 5ft 4 
5 3 ft O 

1 23ft 
OOOO 



. ft 227E * 05 



-POTENTIAL 


X - VELOC I TY 


Y - VELOC I TY 


7ft . OOOO 


- . ft 2 1 4E-02 


- . 1 203E-04 


7 A . OOOO 


- . 82 1 1 E -02 


- . 97A3E- 14 


7ft . OOOO 


- . A202E -02 


. 8 1 1 5E -Oft 


7ft . OOOO 


- . A 1 A AE -02 


0 . 


78 . OOOO 


-.81 ftftE-02 


. 2002E -04 


7ft . OOOO 


- ft 1 4 1 E -02 


- . ft A70E - 1 4 


7 A . OOOO 


- . ft 1 1 2E -02 


. 30ft 1 E -04 


7ft . OOOO 


- . ft077E -02 


0 . 


7ft . OOOO 


- . ft 0 3ft E -02 


. 493 1 E -04 


7ft . OOOO 


- . 7987E-02 


. 1 70 1 E - 1 3 


7ft . OOOO 


- . 7938E -02 


. 7ftft7E -04 


7ft . 0770 


- . A490E - 02 


- . 22ft AE -05 


75 . 0501 


- . ft 4 S4 E - 02 


. 8 305 E -04 


75 . 0894 


- . 84 ftftE-02 


. 1 28 7E -03 


75 . 1019 


- . A447E -02 


. 1 A44E-03 


75.1194 


- . 54 1 AE -02 


24S3E -03 


73.1427 


- . 5353E -02 


32 1 OE -03 


72 . 0570 


- . A7A4E -02 


1 5 1 4E -Oft 


72 . 0579 


- . ft 7 ft 4 E -02 


. 44*AE -04 


72 .0608 


- . 8753E -02 


•456E -04 


72 . 0858 


- . A7A3E -02 


. 1 343E -03 


72 . 0723 


- . A7ft2E -02 


1 854E -03 


72.0783 


- . A7A3E -02 


. 2 Oft 4 E -03 


72 . Oftftft 


- . A7ft4E -02 


2872E -03 


72 . 0935 


- . 5 7 ft ft E -02 


2997E-03 



TOTAL VELO 

. 52 1 4E -02 
5 2 1 1 E *02 

. 5 202 E *02 
ft 1 A9E-02 
ft 1 ftSE -02 
ft 1 4 1 E - 02 
ft 1 1 2E-02 
8077E -02 
A03AE -02 
7957E-02 
7934E -02 
ft 4 ft OE - 02 
84SSE -02 

. 84 BSE -02 
844SE -02 
5422E -02 
83BSE -02 
ft 7 ft 4 E -02 
ft 7 ft 4 E - 02 
87ft4E -02 
ft 7 ft 4 E - O 2 
8754E - 02 
5755E -02 
87ft SE -02 
8790E -02 



2 • 

30 

31 

32 

33 

34 

34 

34 

37 

34 

34 

40 

4 1 

42 

43 

44 

44 

44 

47 

44 

44 

SO 

S 1 

5 2 

63 

44 

S 5 

S 4 

S 7 

54 

44 

60 

6 1 

62 

63 

64 

66 

66 

67 

66 

66 

70 

71 

72 

73 

74 

76 

74 

77 

74 

74 

40 

6 1 

42 

63 

44 

66 

66 

67 

46 

66 

60 

6 1 

62 

63 

64 

66 

6 6 

47 

16 

66 

100 

101 

102 

103 

104 

10S 

106 

1 07 

106 

106 

1 10 

1 1 1 

1 12 

1 1 3 

1 1 4 

1 1 S 

1 1 6 

1 1 7 

1 1 6 

1 16 

1 20 

1 2 1 

1 22 

123 

1 24 

126 

126 

127 

125 

126 

1 30 

1 3 1 

1 32 

133 

134 

135 

136 

1 37 

136 

1 39 

1 40 

1 4 1 

1 42 

143 

144 

146 

146 

147 

1 46 

1 41 

1 60 

1 5 1 

152 

153 

154 

1 55 

156 

157 



140 OOOO 


0 OOOO 


69 3639 


- . 9 1 66E - 02 


-.19326-06 


91666-02 


1 40 . OOOO 


1 7 . OOOO 


693662 


• 9 1 94E -02 


. 10106-03 


. 9 1 856-02 


1 40 . OOOO 


34 . OOOO 


69.3512 


- . 9 1 60E-02 


. 20236-03 


6 1626-02 


1 40 . OOOO 


45 .7836 


69 . 3963 


- . 9 1 44E-02 


. 27196 -03 


9 1 456-02 


1 40 . OOCO 


67 . 8573 


69.4136 


- . 9 1366-02 


. 34136-03 


. 9 1426-02 


140 . OOOO 


69 . 3606 


66 . 4363 


• . 9 1286-02 


. 41276-03 


. 9 1346-02 


• 60.0000 


0 . OOOO 


66 . 6636 


- .95356-02 


. 50366-05 


. 98366-02 


-60.0000 


6 . 6000 


66.5648 


- .66326-02 


. 54716-04 


. 95336-02 


• 60.0000 


1 7 . OOOO 


66.6664 


- . 9 6296 -02 


. 1 1966-03 


. 65306-02 


-60 . OOOO 


28 . 5000 


66 . 6742 


- . 95226 -02 


. 1 6646-03 


. 68236-02 


-60.0000 


34 . OOOO 


68 .6526 


- . 95 1 4e -02 


. 23756-03 


.66 176-02 


-60.0000 


38 . 4367 


66.5691 


- . 96076 -02 


. 26016-03 


. 681 0E-02 


•60.0000 


44 . 5793 


66 . 6967 


- .65006-02 


.31196-03 


. 96056-02 


•60.0000 


60.3190 


66 . 6062 


- . 94926 -02 


. 32606-03 


. 94976-02 


-60.0000 


65 . 7557 


66.6146 


- . 94526 -02 


. 35206-03 


64606-02 


- 60 . OOOO 


61.1663 


66 . 626 1 


- . 94766 -02 


. 40066-03 


. 94846-02 


-SO . OOOO 


66 . 6360 


66 . 6364 


- .94656-02 


. 45276-03 


64796-02 


- 40 . OOOO 


0 . OOOO 


64.6314 


- .97636-02 


- . 698 16 -05 


9763e-02 


-40 . OOOO 


1 7 . OOOO 


64 . 6387 


- . 677ie -02 


. 60966-04 


97766-02 


-40 . OOOO 


34 . OOOO 


64 . 6492 


-.961 66-02 


. 16666-03 


. 96 1 6E -02 


- 40 . OOOO 


44 . 2206 


64 . 6628 


- . 66356 -02 


. 27046-03 


96426-02 


•40 . OOOO 


54.4411 


64 . 6606 


• . 66566-02 


.36646-03 


96636-02 


-40 . OOOO 


64.6617 


64 . 7032 


- . 95566-02 


. 47636-03 


. 95676-02 


0 . OOOO 


0 . OOOO 


62 . 6587 


- . 1 0046-01 


. 36706-04 


. 1 0046-0 1 


0 . OOOO 


6 . 8000 


62 . 6676 


- . 1 0066-01 


- . 22126-04 


. 10066-01 


0 . OOOO 


1 7 . OOOO 


62 . 6568 


- . i oi i e -o i 


.69146-05 


. 101 16-01 


0 . OOOO 


26 . SOOO 


62 . 6546 


- . 10166-01 


- . 54646 -06 


. 101 66-01 


0 . OOOO 


34 . OOOO 


62.6861 


- . 1 023e -01 


. 1 720e -03 


. 10236-01 


0 . OOOO 


36 . 7606 


62.6566 


- . 1 0226-01 


. 74716-04 


. 10226-0 1 


0 . OOOO 


43.6615 


62 . 6567 


- . 10186-01 


. 36036 -03 


. 10196-01 


0 . OOOO 


46 . 3427 


62 . 6625 


- . i oi i e -o i 


. 18666 -03 


.1011 6-01 


0 . OOOO 


53.1236 


62 . 6673 


- . 1 0006-01 


. 50086-03 


. 10026-0 1 


0 . OOOO 


67.6046 


62 . 6760 


- . 69346-02 


.36436-03 


. 99406-02 


O . OOOO 


62 . 6854 


62 . 6647 


- . 95466-02 


. 64656-03 


. 656 66-02 


2 S . OOOO 


0 . OOOO 


61.4413 


- . 93 176-02 


- . 1 7636 -04 


. 9 3 176-02 


26 . OOOO 


1 7 . OOOO 


6 1 .4166 


- .95606-02 


• . 60046 -03 


.95686-02 


25 . OOOO 


34 . OOOO 


6 1 .3636 


- . 1 0326-0 1 


• . 96836 -03 


. 10376-0 1 


26 . OOOO 


43 .0601 


6 1 .3363 


-.10966-01 


- . 67036-03 


. 1 1026-0 1 


26 . OOOO 


62 . 1203 


6 1 . 3276 


-.11716-01 


- . 43256-03 


. 1 1726-01 


26 . OOOO 


61.1 504 


61.3466 


- . 1 2006-0 1 


. 54246-03 


. 1 2026-01 


SO . OOOO 


0 . OOOO 


60 . 3264 


- . 62466-02 


- . 35236-04 


.62466-02 


SO . OOOO 


5 . 6000 


80.3164 


- .53666-02 


- . 70766 -03 


. 63656-02 


SO . OOOO 


1 7 . OOOO 


80 . 2663 


- . 66826-02 


-.13676-02 


. 66766-02 


60.0000 


26 . 6000 


60.1917 


- . 60866-02 


- . 22736-02 


.93766-02 


SO . OOOO 


34 .0000 


60 .0761 


- . 88576-02 


- . 30176-02 


. 10316-01 


60 . OOOO 


36 .2762 


60.0028 


- . 10606 -01 


- . 36706-02 


. 1 1226-01 


50.0000 


42 . 5565 


56.6160 


- . 1 1436-01 


- . 366 1 e-02 


. 1 207E-01 


60 . OOOO 


46 .6377 


56.6286 


- . 12656-01 


- . 40676-02 


. 136 16-01 


50.0000 


61.1166 


56 . 7437 


- . 14606-01 


- . 33336*02 


. 16 176-01 


50.0000 


68 . 3882 


56 . 8624 


- . 16766-01 


- . 18866-02 


. 16876-0 1 


60 . OOOO 


56 . 6764 


56 . 8784 


- . 1 2666 -0 1 


.44566-02 


. 13636-0 1 


64 . 3760 


0 . OOOO 


60 . 1588 


- . 7708e -02 


. 1 6476-04 


. 77056-02 


54 . 3760 


1 7 . OOOO 


60 . 0606 


- . 7841 e -02 


- . 1 580E -02 


. 606 1 6-02 


64 . 3760 


34 . OOOO 


66 . 6740 


- . 66626 -02 


- . 38176-02 


. 68006-02 


54.4167 


42 . 2762 


56 . 6668 


- . 10176-01 


- . 51666-02 


.11416-01 


64 .4663 


60 . 6666 


66.4148 


- . 14476-01 


• .63466-02 


. 18716-01 


S 4 .6000 


66 . 6377 


66 . 1646 


- . 37326-01 


- . 46066-03 


. 3732e-0 1. 


66 . 7600 


0 . OOOO 


56.6623 


- . 73036-02 


. 83026 -06 


. 73036-02 


66 . 7SOO 


8 . 8000 


86.6751 


- . 73766-02 


• . 62666-03 


. 74246-02 


S 6 . 7600 


1 7 . OOOO 


66.6216 


- .74756-02 


-.16616-02 


. 76566*02 


S 6 . 7SOO 


25 . 6000 


56 .8280 


- . 7776e -02 


- . 27636-02 


52526-02 


66.7600 


34 . OOOO 


58 . 6870 


-.61636-02 


- . 36686-02 


90666-02 


S6 . 76 1 7 


36 . OOOO 


56.5880 


- .84666-02 


- . 4676E -02 


. 88386-02 


66 . 6333 


42 . OOOO 


56 . 4645 


- . 88236-02 


- . 58206-02 


. 10836-0 1 


66 . 6760 


46 . OOOO 


56 . 3433 


- . 68766-02 


- . 63856-02 


. 1 2606 0 1 


66.6167 


50 . OOOO 


56.1450 


- . 10 14E -01 


- . 66276-02 


. 13446-0 1 


66 . 6663 


64 . OOOO 


56 . 6306 


- . 17606-01 


- . 26446 -01 


. 33456-01 


56 . OOOO 


6 6 . OOOO 


66 . OOOO 


- . 3 1 466 -01 


- . 44426 -01 


54436-01 


80. 3750 


0 . OOOO 


56 . 9336 


-.71 60e -02 


. 1 4736-04 


. 7 1 606-02 


60 . 3760 


1 7 . OOOO 


56.6617 


-.731 4E -02 


• . 1 6766 -02 


76046-02 


80 . 3760 


34 . OOOO 


56 . 6220 


- . 76486-02 


- . 40326 -02 


. 5621 6-02 


80 . 2800 


42 . OOOO 


96 . 4240 


- . 63236-02 


• . 66486-02 


. 10236-01 


60 . 1260 


50 . OOOO 


66 .0806 


- . 68306-02 


-.51 186-02 


. 1 2076-01 


80.0000 


66 . OOOO 


56 . OOOO 


• . 1 6626 - 1 2 


• . 46626 -01 


. 46626-0 1 


82 . OOOO 


0 . OOOO 


58 . 6760 


- . 68626-02 


. 1 8366-04 


. 86826-02 


82 . OOOO 


6 . 6000 


56 . 6561 


- . 70456-02 


• .86636-03 


. 70676-02 


62 . OOOO 


1 7 . OOOO 


56 . 6030 


- . 71226-02 


- . 1 7046-02 


. 73236-02 


62 . OOOO 


25 . 6000 


68 . 7067 


• . 731 16-02 


• . 26646 -02 


. 76526-02 


62 . OOOO 


34 . OOOO 


58 . 5666 


- . 76276-02 


- . 40666-02 


. 65686-02 


61 .8333 


36 . OOOO 


66 . 4736 


- .76626-02 


- . 6 1446-02 


.92266-02 


61 . 6667 


42 .0000 


56 . 3688 


- . 78416-02 


- . 60166-02 


. 88836 -02 


61 . SOOO 


48 .0000 


56 . 2308 


- . 7581 e -02 


- . 66826-02 


. 1 1446-0 1 


61 . 3333 


50 . OOOO 


66 . 0371 


- .68656-02 


- . 66636-02 


. 1 2286-0 1 


61 . 1687 


84 . OOOO 


66 . 7428 


.926 IE -03 


- . 26686-01 


. 258 16-0 1 


8 1 . OOOO 


86 . OOOO 


56 . OOOO 


. 22636 -01 


- . 42186-01 


. 47876-01 


64 . OOOO 


0.0000 


56 . 6070 


- . 86026-02 


. 13646-04 


. 66026-02 


64 . OOOO 


1 7 . OOOO 


56 . 7328 


- . 86026-02 


- . 17336-02 


.711 66-02 


64 . OOOO 


34 . OOOO 


86.4662 


- . 71426-02 


- . 4 1426-02 


. 62566-02 


83 . 6766 


42.1010 


59 . 2663 


- . 70846-02 


- . 60666-02 


. 63246-02 


63.3512 


50 . 202 1 


66 . 8732 


- . 44626-02 


- . 86246-02 


. 10596-01 


63.0288 


66 . 3031 


58 . 4487 


. 31466 -OI 


- . 1 050E -OI 


. 33 1 6 6-0 1 


86 . OOOO 


0 . OOOO 


56 . 7398 


- . 88886-02 


. 1 4146*04 


. 669 66-02 


66 . OOOO 


6 . 6000 


68.7216 


- . 6 82 6 E - 02 


• . 55156-03 


. 66646-02 


66 . OOOO 


1 7 . OOOO 


58.6650 


- . 88736-02 


• . 1 765E -02 


. 69006-02 


66 . OOOO 


26 . 6000 


59.5883 


- . 67366-02 


• . 262 1 6-02 


. 73426-02 


66 . OOOO 


34 . OOOO 


59.4167 


- . 87796-02 
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1 7 . OOOO 


55.4516 


- . 4540E-03 


- . 794 6E -03 


. 9 1 52E-03 


463 


146. OOOO 


26 . 5000 


58.4102 


- . 4 1 37E -03 


-.1181 E -02 


. 1 232E-02 


464 


146. OOOO 


34 . OOOO 


58 . 3533 


- . 3547E-03 


-.151 3E -02 


. 1 557E -02 


465 


146. 1560 


36 .0265 


58.3214 


- . 3345E -03 


- . 1 824E -02 


. 1 658E -02 


466 


146.3100 


42 . 059 1 


55 . 2673 


- . 3086E -03 


- . 1 7 4 4 E * 02 


. 1 77 IE -02 


467 


146 . 4646 


46 . 0696 


55.2519 


- . 3470E -03 


- . 1 7 1 6E -02 


. 1 753E -02 


466 


146.6169 


50 . 1191 


55.2176 


- . 4464E -03 


- . 1 665E-02 


. 1 743E -02 


466 


1 46. 7749 


54 . 1 477 


56 . 1863 


- . 7386E -03 


- . 9733E-03 


. 1 222E -02 


470 


146.6266 


59 . 1 772 


56 . 1772 


- . 671 6 E - 03 


- . 5391 E-04 


. 6752E -03 


471 


145. OOOO 


0 . OOOO 


55.4612 


- . 3 6 6 6 E - 03 


. 1005E-OS 


. 3666E -03 


472 


146 . OOOO 


1 7 . OOOO 


55 . 4476 


- . 3600E -03 


- . 7902E-03 


. 6653E-03 


473 


145 . OOOO 


34 . OOOO 


56 . 3500 


- . 2632E -03 


- . 1 50 1 E -02 


. 1526E-02 


474 


146 . 3067 


42 . 0530 


56 . 2644 


- . 2732E-03 


- . 1 7 4 6 E - 02 


. 1768E-02 


475 


145.6195 


50 . 1061 


56.2125 


- . 5 1 72E -03 


- . 1 674E-02 


. 1 6 4 4 E - 02 


476 


146.6262 


56.1561 


56.1641 


-.161 3 E - 02 


- . 2425E-03 


. 1 626E -02 


477 


1 50 . OOOO 


O . OOOO 


56 . 4780 


- . 2765E -03 


- . 7736E-07 


. 2765E-03 


476 


160. OOOO 


6 . 6000 


56 . 4667 


- . 2744E-03 


-.361 IE-03 


. 4778E-03 


476 


160. OOOO 


1 7 . OOOO 


56 . 4447 


• . 2700E -03 


- . 7666E-03 


. 63 1 6E-03 


460 


160 OOOO 


26 . 5000 


56 . 4036 


- . 2475E -03 


-.11 45E-02 


.117 IE-02 


46 1 


160. OOOO 


34 . OOOO 


56 . 3474 


- . 2210E-03 


- . 1 462E-02 


. 1 806E -02 


462 


160. 1646 


36 . 0236 


56.3160 


- . 21 6 1 E -03 


-.16 1 6E -02 


. 1 633E -02 


463 


1 50 . 3066 


42 . 0470 


56.2620 


- . 2122E-03 


- . 1 76 6E -02 


. 1 766E -02 


464 


1 50 . 4643 


46 . 0705 


56 . 2454 


- . 3546E -03 


- . 1 6 3 6 E * 0 2 


. 1 873E -02 


466 


150.6160 


60 . 0640 


56 . 2074 


- . 4676E-03 


- . 2066E -02 


. 21 46E -02 


466 


1 50 . 7739 


54.1176 


56.1661 


- . 1 6 1 7E-02 


- . 1 582E -02 


. 2406E -02 


467 


150.6265 


56 . 14 10 


56.1410 


- . 276 2E-02 


- . 3342E -03 


. 26 1 2E -02 


466 


152. OOOO 


0 . OOOO 


68 . 4757 


- . 1 637E -03 


- . 4636E-06 


. 1 637E -03 


466 


152. OOOO 


1 7 . OOOO 


56 . 4426 


- . 1 760E-03 


- . 7640E -03 


. 604 1 E -03 


460 


152 . OOOO 


34 . OOOO 


56 . 3456 


- . 1 462E-03 


-.1 46 6E -02 


. 1463E -02 



46 1 


162.3214 


42 .0235 


58 . 2802 




- . 1 630E -03 


- . 1 757E -02 


1764E-02 


462 


1 52 . 6426 


60 . 0470 


59.2013 




- . 6024E -03 


- . 206 1 E -02 


21 66E-02 


463 


1 62 .6643 


56 .0706 


59 . 0677 




- . 7006E -02 


- . 2307E -02 


7376E -02 


464 


1 54 . OOOO 


O . OOOO 


59 . 4743 




- . 6226E -04 


- . 6246E -06 


6226E -04 


466 


1 54 . OOOO 


8 .6000 


56 .4660 




-.61 26E -04 


- . 3600E -03 


4006E -03 


46 6 


1 54 . OOOO 


1 7 . OOOO 


56.4412 




- . 6650E-04 


- . 7624E -03 


7876E-03 


467 


154 . OOOO 


25 . 5000 


56 . 4003 




- . 6240E-04 


- . 1 1 35E -02 


1 141E-02 


466 


154 . OOOO 


34 .OOOO 


58 .3444 




- . 741 2E-04 


- . 1 48 1 E-02 


1 483E -02 


466 


154.1667 


36 .OOOO 


56.3134 




- . 71 20E-04 


- . 1 630E-02 


1 632E -02 


5 00 


154 . 3333 


42 .OOOO 


66 . 279 1 




- . 6423E -04 


- . 1 746E -02 


1 7SOE -02 


601 


1 54 . 5000 


46 .OOOO 


58.2415 




- . 1323E-03 


- . 2067E -02 


207 IE-02 


502 


154.6667 


50.0000 


58.1662 




. 3066E -04 


- . 1 782E-02 


1 763E-02 


503 


154 .6333 


54 . OOOO 


56 . 1367 




- . 1 446E -02 


- . 4645E-02 


5 05 7 E - 02 


504 


155 . OOOO 


56 . OOOO 


58 . OOOO 




- . 4660E-02 


- . 76 6 5 E -02 


6372E-02 


505 


1 56 . OOOO 


O . OOOO 


56 . 4740 




- . 4662E-04 


- . 1044E-05 


4594E -04 


506 


155. OOOO 


1 7 . OOOO 


56 . 4406 




- . 4453E-04 


- . 7622E-03 


7635E-03 


507 


165. OOOO 


34 . OOOO 


58 . 3442 




- . 3623E-04 


- . 1 482E-02 


1482E-02 


506 


155.1 667 


42 . OOOO 


58 . 2780 




- . 1 655E-04 


- . 1756E-02 


1 7 5 6 E - 02 


506 


1 55 .3333 


50 . OOOO 


56 . 16 67 




. 98 2 1 E - 04 


- . 1728E-02 


1 730E-02 


5 1 O 


155 . 5000 


56 . OOOO 


56 . OOOO 




. 2726E- 1 2 


- 576SE-02 


8 7 8 6 E • 02 


6 1 1 


156. OOOO 


0 . OOOO 


56 . 4739 




- . 1 306E-06 


- . 1 1 1 IE -05 


1 1 1 8E-05 


5 1 2 


156. OOOO 


8 . 5000 


55 . 4656 




. 1046 E-OS 


- . 3666E-03 


3869E -03 


5 1 3 


156 . OOOO 


1 7 .OOOO 


56 . 4407 




- . 4220E-06 


- . 7621E-03 


7521E -03 


5 1 4 


156. OOOO 


26 . 6000 


55 . 3669 




. 1669E-OS 


- . 1 1 37E -02 


1 1 37E -02 


5 1 5 


166 . OOOO 


34 . OOOO 


56 .344 1 




. 30S5E-07 


- . 1483E-02 


1 463E -02 


5 1 6 


156. OOOO 


38 . OOOO 


56.3130 




. 1 766E-05 


-.1S2SE-02 


1 625E-02 


5 1 7 


156. OOOO 


42 . OOOO 


55 . 2790 




. 2665E -05 


- . 1 765E -02 


1 78SE-02 


5 1 6 


156. OOOO 


46 . OOOO 


58 .2409 




. 201 6E -04 


- . 2052E-02 


2053E -02 


5 1 9 


156. OOOO 


50.0000 


55.1969 




. 6436E -04 


- . 1 775E -02 


1 750E -02 


520 


156. OOOO 


54 . OOOO 


65 . 1358 




. 1 1 97E -03 


- . 4 9 2 2 E - 02 


4923E -02 


52 1 


156. OOOO 


66 . OOOO 


56 . OOOO 




-.31 58E -04 


- . 5656E-02 


5659E -02 


SEEPAGE PACE CDMPU TAT 1 


ION RESULTS 












N 5 * 


66 HE « 113 


AVE. VEL > 


. 51200E-01 


D 6 


> . 20000E +0 1 


DO • . 10240E+00 






TDTAL DISCHARGE 


> . 1 0 2 4 OE * 00 










NB ■ 


113 NE « 196 


AVE. VEL * 


.361 07E -01 


D5 


* . 2 2000E +0 2 


DO * . 5 3 6 3 6 E + 00 






TDTAL DISCHARGE 


« . 6 407 6E +00 










NB > 


195 NE > 215 


AVE. VEL = 


. 26362E-01 


D5 


t . 20000E ♦ 0 1 


DO « . S2704E -Ol 






TDTAL DISCHARGE 


> . 6 6 3 4 6 E +00 










N 6 « 


215 NE > 300 


AVE. VEL > 


. 20451 E - O 1 


D5 


s , 2 2000 E +02 


DO > . 4 5 05 8 E +00 






TDTAL DISCHARGE 


. 1 4 4 4 OE * O 1 










N 5 > 


300 NE ■ 317 


AVE. VEL « 


. 15646E-01 


DS 


« . 20000 E +0 1 


DO « . 3 1 29 1 E - O 1 






TDTAL DISCHARGE 


« . 1 4753E + 01 










N 6 * 


317 NE > 402 


AVE. VEL > 


. 1 2764E-01 


D5 


• . 2 2000E +0 2 


DO > .2812SE+00 






TDTAL DISCHARGE 


« . 1 75 6 6 E +0 1 










N6 > 


402 NE ■ 419 


AVE. VEL > 


. 10683E -01 


DS 


> . 2 OOOOE +01 


DO > .21 1 66E-01 






TDTAL DISCHARGE 


s . 1 777 6 E +0 1 










N 6 « 


419 NE = 504 


AVE. VEL > 


. 985 10E-02 


DS 


= . 2 2000E +0 2 


DO b . 2 1 7 35 E *00 






TDTAL DISCHARGE 


« . 1S8S1E+01 










N 6 = 


604 NE * 621 


AVE. VEL ■ 


. 901 57E-02 


D 6 


• . 1 OOOOE +01 


DO « . 80 1 5 7E - 02 






TDTAL DISCHARGE 


= . 2004 2E +0 1 











■ SP • 25 . . K • .2..D • 20 FT. ..END OF PROS LEM. 
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